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Abstract. These are extended notes of a series of lectures given at tfie Xlth Brazilian School of Cosmology and Gravitation. 
They provide a selection of topics at the intersection of particle and astrophysics. The first part gives a short introduction 
to the theory of electroweak interactions, with specific emphasize on neutrinos. In the second part we apply this framework 
to selected topics in astrophysics and cosmology, namely neutrino oscillations, neutrino hot dark dark matter, and big bang 
nucleosynthesis. The last part is devoted to ultra high energy cosmic rays and neutrinos where again particle physics aspects 
are emphasized. The often complementary role of laboratory experiments is also discussed in several examples. 
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1. INTRODUCTION AND REMINDER: FERMI THEORY OF WEAK INTERACTIONS 



Good introductory texts on particle physics are contained in Ref. [1] (more phenomenologically and experimentally 
oriented) and in Refs. [2, 3, 4]. Here we will only recall the most essential facts. 

We will usually use natural units in which h ~ c ^ k = 1 , unless these constants are explicitly given. 
Neutrinos only have weak interactions. Historically, experiments with neutrinos obtained from decaying pions and 
kaons have shown that charged and neutral leptons appear in three doublets: 

TABLE 1. The lepton doublets 
q Le = l iju = 1 if = 1 



Charge q and lepton numbers Le, L^, and are conserved separately. There are corresponding doublets of anti- 
leptons with opposite charge and lepton numbers, denoted by v,- for the anti-neutrinos and by the respective positively 
charged anti-leptons. 

Therefore, allowed reactions include n pe^Ve (nuclear j3 -decay), VeP ne^ (inverse neutron decay) ;r+ ^ 
M^V|i, V^n pjj.^, but exclude Vn/? nji^, /i+ — > e+y. 

The "neutrino" is thus defined as the neutral particle emitted together with positrons in j3+-decay or following 
K-capture of electrons. The "anti-neutrino" accompanies negative electrons in /3~ -decay. 

Lifetimes for weak decays are long compared to Ufetimes associated with electromagnetic (~ 10" '^s) and strong 
(~ 10"^^ s) interactions. A weak interaction cross section at ~ IGeV interaction energy is typically ~ lO''^ times 
smaller than a strong interaction cross section. 

Weak interactions are classified into leptonic, semi-leptonic, and non-leptonic interactions. 

Fermi's golden rule yields for the rate F of a reaction from an initial state i to a final state / the expression 

where Mij = (/|//int) is the matrix element between initial and final states i and / with Hi^i the interaction energy, and 
dN/dEf is the final state number density evaluated at the conserved total energy of the final states. 
As an example, we compute the rate F for inverse /3 -decay 

Ve + p^n + e+. (2) 

We use the historical Fermi theory after which such interactions are described by point-like coupUngs of four fermions, 
symbolically = Gp / d^xy/^, with Fermi's coupling constant Gp. This yields 

2k T: .T dN 

F=— G^Mp— , (3) 
n dEf 

where symbolically M — f d^Xj/^ which incorporates the detailed structure of the interaction. If we normalize the 
volume V to one, M is dimensionless and of order unity, otherwise M scales as V"'. In fact, it is roughly the spin 
multiphcity factor, such that |Mp ~ 1 if the total leptonic angular momentum is 0, thus involving no change of spin in 
the nuclei ("Fermi transitions"), whereas |Mp ~ 3 if the total leptonic angular momentum is 1, thus involving a change 
of spin in the nuclei ("Gamow-Teller transitions"). The final state density of a free particle is 

(4) 



Therefore, taking into account energy-momentum conservation, we get in the center of mass (CM) frame the phase 
space factor for the two body final state 

.A. _ f y.3p, y,3p^ (^-m^SiE. + E.-Eo), (5) 



dEf J {2nhf {2nh)^ V 

where pe, p„, £„, are momenta and kinetic energies of the electron and the final state nucleus, respectively, and 
Eq is the total initial energy. Integrating out one of the momenta gives Pf = Pe = Pn so that energy conservation 
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Eo = {Pf + {p^f + m^yl^ gives the factor dpf/dEo = {pf/Eg + Pf/En) ^ = ^ with v/ being the relative 
velocity of the two final state particles. This yields 

^ = JLp1 (6) 

dEf 2n^ Vf' ^ ' 

We are now interested in the cross section a of the two-body reaction Eq. (2) defined by 

r = anm , (7) 

where n,- = and v, are density and velocity, respectively, of one of the incoming particles in the frame where the 
other one is at rest. Putting this together with Eqs. (3) and (6) finally yields 

a{VeP^ne+) = ^\Mif\^-^. (8) 

K ViVf 

For Pf—l MeV this cross section is ~ 10~^^ cm^. In a target of proton density Hp this gives a mean free path defined 
by lvnpa{Vep nc+) ^ 1. For water this turns out to be ~ 30pc which demonstrates the experimental challenge 
associated with detection of neutrinos. 

The first detections of this reaction was made by Reines and Cowan in 1959. The source were neutron rich fission 
products undergoing /3 -decay n ^ pe~Ve-A 1000 MW reactor gives a flux of ^ lO'^cm^^s^' V^s which they observed 
with a target of CdCl2 and water. Observed are fast electrons Compton scattered by annihilation photons from the 
positrons within ^ 10"^ s of the reaction ("prompt pulse") y— rays from the neutrons captured by the cadnnium about 
10"^ s after the reaction ("delayed pulse"). 



2. DIRAC FERMIONS AND THE V-A INTERACTION 

Given the experimentally established fact that electroweak interactions only involve left-handed neutrinos we now 
want to work out the detailed structure of these interactions. In order to do that we first have to introduce the Dirac 
fermion. 



2.1. Dirac Fermions as Representations of Space-Time Symmetries 

The Poincare group is the symmetry group of special relativity and consists of all transformations leaving invariant 
the metric 

ds^ = -{dxY + {dx^f + {dxY + [dx^f , (9) 
where is a time coordinate and x^,x^, and x^ are Cartesian space coordinates. These transformations are of the form 

x'^' = A'ix''+a^', (10) 

where defines arbitrary space-time translations, and the constant matrix Ay satisfies 

T7^vA|^A^ = T7pCT , (11) 

where rjjxv = diag(— 1,1,1,1). The unitary transformations on fields and physical states induced by Eq. (10) satisfy 
the composition rule 

J/(A2,a2)f/(Ai,ai) =J/(A2Ai,A2fli+fl2). (12) 

Important subgroups are defined by all elements with A = 1 (the commutative group of translations), and by all 
eleraentswiiha'^ = 0[thehomogeneousLorentzgroup S0{3, 1) of matrices Ay satisfying Eq. (11)]. The latter contains 
the subgroup 50(3) of all rotations for which Aq = 1, A^ = A^ = for /j = 1 , 2, 3. 

The general infinitesimal transformations of this type are characterized by an anti-symmetric tensor £»{f and a vector 

A(J = 5i' + a)if a^' = e^'. (13) 
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Any element f/(l + co, e) of the Poincare group which is infinitesimally close to the unit operator can then be expanded 
into the corresponding hermitian generators J^^ and P^, 

U{l + (0,e) = l + ^i(0^,vJ'''' -ie^P^ . (14) 
It can be shown that these generators satisfy the commutation relations 

i[P^,JP''] = j^MPp^ _ j^M^^pP (15) 

[P^',P''] = 0. 

The P^ represent the energy-momentum vector, and since the Hamiltonian H = P^ cormnutes with the spalial pseudo- 
three-vector J = {J^^,J^^,J^'^), the latter represents the angular-momentum which generates the group of rotations 
50(3). 

The homogeneous Lorentz group implies that the dispersion relation of free particles is of the form 

E^{p)^p^+M^. (16) 

for a particle of mass M, momentum p, and energy E. If one now expands a free charged quantum field \^(x) into its 
energy-momentum eigenfunctions and interprets the coefficients a (p) of the positive energy solutions as annihilator of 
a particle in mode p, then the coefficients ^^(p) of the negative energy contributions have to be interpreted as creators 
of anti-particles of opposite charge, 

Va(x)= £ a(p)M(p)e-'^W'+'P-''+ £ ^^(p)v(p)e'^(P)'-*-''. (17) 

p,E(p)>0 P,£(p)<0 

Canonical quantization, shows that the creators and annihilators indeed satisfy the relations, 

fl,-(p),4(P')l - k-(p),^J(p')l =5,v5(p-p'), (18) 

L J ± L J ± 

where i, i' now denote internal degrees of freedom such as spin, and [., .]± denotes the commutator for bosons, and the 
anti-commutator for fermions, respectively. 

Fields and physical states can thus be characterized by their energy-momentum and spin, which characterize their 
transformation properties under the group of translations and under the rotation group, respectively. Let us first focus 
on fields and states with non-vanishing mass. In this case one can perform a Lorentz boost into the rest frame where 
P^ = (M, 0,0,0) with M the mass of the state. P^ is then invariant under the rotation group 50(3). The irreducible 
unitary representations of this group are characterized by a integer- or half-integer valued spin j such that the 27 + 1 
states are characterized by the eigenvalues of 7, which run over — 7 , — 7 ' + 1 , • • • , 7 — 1 , y . Note that an eigenstate with 
eigenvalue a of is multiplied by a phase factor g^^itr under a rotation around the «— axis by 2n, and a half-integer 
spin state thus changes sign. Given the fact that a rotation by 2k is the identity this may at first seem surprising. 
Note, however, that normahzed states in quantum mechanics are only defined up to phase factors and thus a general 
unitary projective representation of a syrmnetry group on the Hilbert space of states can in general include phase 
factors in the composition rules such as Eq. (12). This is indeed the case for the rotation group 50(3) which is 
isomorphic to ^3 /Z2, the three-dimensional sphere in Euclidean four-dimensional space with opposite points identified, 
and is thus doubly connected. This means that closed curves winding n times over a closed path are continuously 
contractible to a point if n is even, but are not otherwise. Half-integer spins then correspond to representations for 
which U(A\ )U(A2) = (— )"C/(Ai A2), where n is the winding number along the path from 1 to Ai, to A1A2 and back 
to 1, whereas integer spins do not produce a phase factor. 

With respect to homogeneous Lorentz transformations, there are then two groups of representations. The first one 
is formed by the tensor representations which transform just as products of vectors, 

Wy":: =A''pA^---w§:::. m 

These represent bosonic degrees of freedom with maximal integer spin 7 given by the number of indices. The simplest 
case is a complex spin-zero scalar (j) of mass m whose standard free Lagrangian 

^^ = ~{^^,(^^^^'(^-n^(^^(^), (20) 
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leads to an equation of motion known as Klein-Gordon equation, 

{d^d^ -n?)^=0. (21) 
In the static case = Q this leads to an interaction potential 

-mr 

y{r)=g\g2 , (22) 

r 

between two "charges" gi and g2 which correspond to sources on the right hand side of Eq. (21). The potential for 
the exchange of bosons of non-zero spin involve some additional factors for the tensor structure. Note that the range 
of the potential is given by ~ mT^ . In the general case ^0 the Fourier transform of Eq. (21) with a delta-function 
source term on the right hand side is —i/{p^ + m^). A four-fermion point-like interaction of the form GpXj/^ can 
thus be interpreted as the low-energy limit p^ ^ of the exchange of a boson of mass m. Later we will realize that 
the modern theory of electroweak interactions is indeed based on the exchange of heavy charged and neutral "gauge 
bosons". In the absence of sources, Eq. (21) gives the usual dispersion relation =p^ -\- for a free particle. 

The second type of representation of the homogeneous Lorentz group can be constructed from any set of Dirac 
matrices satisfying the anti-commutation relations 

{ty] = 2^^'\ (23) 
also known as Clifford algebra. One can then show that the matrices 

indeed obey the commutation relations in Eq. (15). The objects on which these matrices act are called Dirac spinors 
and have spin 1 /2. In 3-h1 dimensions, the smallest representation has four complex components, and thus the are 
4x4 matrices. A possible representation of Eq. (23) is 

ri=(Z i=h2,3, ro = *fn (25) 



id ' ' ' \o -1 

where cr, are the Pauli matrices. 

The standard free Lagrangian for a spin- 1/2 Dirac spinor Xj/ of mass m, 

^v = -V{fd^+m)\i/, (26) 

where = V^^^V . leads to an equation of motion known as Dirac equation, 

(/a^+m)vA = 0. (27) 

Its free solutions also satisfy the Klein-Gordon equation Eq. (21) and are of the form Eq. (17) where, up to a 
normaUzation factor A'^, 



u{p)=N{ ^ , V(p) =N[E+m'\. (28) 



Here, u and v are 4-spinors, whereas U and v are two-spinors. 
It is easy to see that the matrix 



Y5 = -ifYV)^ = -(l J) 



(29) 

is a pseudo-scalar because the spatial Y change sign under parity transformation, and satisfies 

)^ = i {/,r5} = o [y'^^r5] = o. (30) 

A four-component Dirac spinor \j/ can then be split into two inequivalent Weyl representations xj/l and Xj/R which are 
called left-chiral and right-chiral, 

1 + y«: 1 — y<i 

Y=WL + WR = ^r^W+^r^W- (31) 
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Note that according to Eqs. (30) and (31) the mass term in the Lagrangian Eq. (26) flips chirality, whereas the kinetic 
term conserves chirality. 

The general irreducible representations of the homogeneous Lorentz group are then given by arbitrary direct 
products of spinors and tensors. We note that massless states form representations of the group S0{2) leaving invariant 
P^, instead of of 5(9(3). The group 5(9(2) has only one generator which can be identified with helicity, the projection 
of spin onto three-momentum. For fermions this is the chirahty defined by 75 above. 

In the presence of mass the relation between chirality and helicity H = a p/pis more comphcated: 



i±r5 



m(p) 



l±75 



V(P) 



N 
2 

N 
2 

N 
^2 



It 



It 



g p 

E + m 

P 



E 

It 



u 

1+H 



(32) 



E + m 



1±- 



H 



■ + 



E 
1±- 



1-H 
2 

1-H 



u 

+u 



E + m 



V 

TV 



From this follows that in a chiral state utfi the helicity polarization is given by 



jLJi _ jL,R 



Pl,r = ,0 

/+'+/: 



(33) 



where Z^** are the intensities in the H = ±l states for given chirahty L or R. Note that due to Eq. (17) the physical 
momentum of anti-particles described by the v spinor is — p in this convention, and therefore the helicity polarization 
for anti-particles in pure chiral states are opposite from Eq. (33): Left chiral particles are predominantly left-handed 
and left-chiral anti-particles are predominantly right-handed in the relativistic limit. Furthermore, helicity and chirality 
commute exactly only in the limit p 3> m, v ^ 1. The experimental fact that observed electron and neutrino helicities 
are +v for particles and anti-particles, respectively, where v is the particle velocity, now imphes that both electrons 
and neutrinos and their anti-particles are fully left-chiral. 



2.2. The V-A Coupling 

Since Dirac spinors have 4 independent components, there are 16 independent bihnears listed in Tab. 2. Using the 
equality 

Yl = 7brnro, (34) 

which can easily be derived from Eq. (25), one sees that the phase factors of the bihnears in Tab. 2 are chosen such 
that their hermitian conjugate is the same with v^i <-^- 1/^2- 

TABLE 2. The Dirac bilinears. For 

Vi = V2 these are real. 





scalar 


5 




4-vector 


V 




tensor 


T 




axial 4-vector 


A 




pseudo-scalar 


P 



Lorentz invariance implies that the matrix element of a general j3 -interaction is of the form 

M = Gf £ C,(VAi^,T2)(VA3^,T4), (35) 

such that only the same types of operators from Tab. 2 couple and common Lorentz indices are contracted over. 

Eq. (35) is a Lorentz scalar However, we know that electroweak interactions violate parity and thus we have to add 
pseudo-scalar quantities to Eq. (35). Equivalently, we can substitute any lepton spinor in Eq. (35) by j(l + 75)1/- 
This is correct at least for the interactions with charge exchange, the so called charged current interactions, for which 
we know experimentally that both neutrinos and charged leptons are fully left-chiral. Using Eq. (30), this leads to 

terms of the form 

/i;L^v,-,L = /r(l- 75)^(1 + 75)Vi, / = e,M,T, (36) 



Lectures on Astroparticle Physics 



February 2, 2008 



7 



which implies that only the V and A type interactions from Tab. 2 can contribute. The general form of charged current 
interactions involving neutrinos is therefore usually written as 



(37) 



3. DIVERGENCES IN THE WEAK INTERACTIONS AND RENORMALIZABILITY 

A incoming plane wave yi, = e'*^^ of momentum k in the z-direction can be expanded into incoming and outgoing radial 
modes e"'*'' and e'^'', respectively, in the following way 

e*^ = :^7rL(2/+ 1) [(-l)'e-'*^'"-e'H P;(cos0), (38) 

where Pi{x) are the Legendre polynomials and cos0 =z/r. Scattering modifies the outgoing modes by multiplying 
them with a phase e^'^' and an ampUtude r]i with < tJ; < 1. The scattered outgoing wave thus has the form 

"Kscatt = -^^£(2^ + 1)-^^^^ f/(cos0) = —F{e), (39) 

where is called the scattering amplitude. 

Let us now imagine elastic scattering in the CM frame, where momentum and velocity v are equal before and 
after scattering. The incoming flux is then v| i/Ajp = v and the outgoing flux through a solid angle dSl is v| i/Ascatt '\^r^d€l = 
v\F{Q)'\^dQ.. The definition Eq. (7) of the scattering cross section then yields 



= 1^(0)1 • (40) 

el 



Using orthogonality of the Legendre polynomials, / dClPi {£l)Pi' {CI) = 47r5;// /(2/ + 1 ), in Eq. (39), we obtain for the 
total elastic scattering cross section 



(Te, = ^£(2/ + 1) 

P* I 



2i 

For scattering of waves of angular momentum / this results in the upper limit 



(41) 



4n 

C7eM<— (2/+1), (42) 

p* 

which is called partial wave unitarity. 

On the other hand, in Fermi theory typical cross sections grow with as in Eq. (8) and violate Eq. (42) for s-waves 
(/ = 0) for 

— ^>— , (43) 
n pi 

where we used — 4 for the sum over polarizations. This occurs for > (tt/Gf)'''^ — 500 GeV. Such energies 

are nowadays routinely achieved at accelerators such as in the Tevatron at Fermilab. As will be seen in the next section, 
this is ultimately due to the fact that the coupling constant Gp has negative energy dimension and corresponds to a 
non-renormalizable interaction. This will be cured by spreading the contact interaction with the propagator of a gauge 

-1/2 

boson of mass M ~ Gp ' ~ 300 GeV. This corresponds to multiplying the l.h.s. of Eq. (43) with the square of the 
propagator, ~ (1 + p^/M^)~^, which thus becomes 4/ {npl) for — > oo. This scaling with is of course a simple 

consequence of dimensional analysis. As a result, partial wave unitarity is not violated any more at high energies, at 
least within this rough order of magnitude argument. The gauge theory of electroweak interactions discussed below is 
renormalizable. 

Theories which contain only coupUng terms of non-negative mass dimension lead to only a finite number of graphs 
diverging at large energies. It turns out that these divergences can be absorbed into the finite number of parameters of 
the theory which is why they are caUed renormalizable. 
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Good examples of non-renormalizable interaction terms are given by 

— xj/[Y^„Yy]xifF^\ —xjfY,[Y^,Y,]xifF^\ (44) 

where M is some large mass scale presumably related to grand unification and e is the (positive) electric charge 
unit. The gauge invariant field strength tensor F^v = d^Ay — d^uAf^ in terms of the gauge potential Eq. (57) below 
represents the electric field strength £' = —doA^ — diA^ = F"' and magnetic field strength B' = s'j^djAii = s'-'^Fji^ 
where Latin indices represent spatial indices and is totally anti- symmetric with e^^-' = 1. As a consequence, in 
the non-relativistic limit, Eq. (44) reduce to a magnetic and electric dipole moment of the Xj/ field, respectively, of size 
4e/M. Note that these are even and odd, respectively, under parity and time reversal. 

Other non-renormalizable terms may arise from Lorentz symmetry violation by physics close to the grand unifica- 
tion scale M. In Sect. 9.5 we wiU see how high energy astrophysics can constrain such terms and thus physics beyond 
the Standard Model to precisions greater than laboratory experiments. 



4. GAUGE SYMMETRIES AND INTERACTIONS 
4.1. Symmetries of the Action 

Lorentz invariance suggests that the action should be the space-time integral of a scalar function of the fields y/i{x,t) 
and their space-time derivatives dn\l/i{x,t), and thus that the Lagrangian should be the space-integral of a scalar called 

the Lagrangian density .if, 

S[\i/] = [d*x^[w,{x),df,xi/i{x)], (45) 



where x = (x,f) from now on. In this case, the equations of motion read 

which are called Euler-Lagrange equations and are obviously Lorentz invariant if .if is a scalar. 

Symmetries can be treated in a very transparent way in the Lagrangian formalism. Assume that the action is invari- 
ant, 5S = 0, independent of whether V'i (x) satisfy the field equations or not, under a global symmetry transformation, 

5 Xj/i (x) = ie^i [yj (x) , dn (x)] , (47) 

for which e is independent of x. Here and in the following expUcit factors of i denote the imaginary unit, and not an 
index. Then, for a space-time dependent e(x), the variation must be of the form 

SS = -J d'^xJf' [x, Xj/j (x) , Xj/j W] a^e (jc) . (48) 

But if the fields satisfy their equations of motion, 8S ^0, and thus 

which implies Noethers theorem, the existence of one conserved current for each continuous global symmetry. If 
Eq. (47) leaves the Lagrangian density itself invariant, an explicit formula for J>^ follows immediately, 

J^' = -i^r^^^i, (50) 

where we drop the field arguments from now on. 

As opposed to a global symmetry, Eq. (47), which leaves a theory invariant under a transformation that is the same 
at all space-time points, a gauge symmetry is more powerful as it leaves invariant a theory, i.e. dJf' = 5S = 0, under 
transformations that can be chosen independently at each space-time point. Gauge symmetries are usually also linear 
in the (fermionic) matter fields which we represent here by one big spinor ^{x) that in general contains Lorentz spinor 
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indices as well as some internal group indices on which the gauge transformations act. For real infinitesimal e" (x) we 
write 

5v/(x) = ;e"(x)fav/(x) , (51) 

where a labels the different independent generators ta of the gauge group. A finite gauge transformation would be 
written as ^{x) — > exp(jA"(x)?(;e)vW and reduces to Eq. (51)in theUmit A"(x) = e"(x) — > 0. The hermitian matrices 
ta form a lie algebra with commutation relations 

[ta,t^] = iCl^ty, (52) 
where the real constants C^n are called structure constants of the Lie algebra, and are anti- symmetric in a/3. 



4.2. Gauge Symmetry of Matter Fields 



If the Lagrangian contained no field derivatives, but only terms of the form i/a(x)^ • • • V'^(x), there would be no 
difference between global and local gauge invariance. However, dynamical theories contain space-time derivatives 
d^'^f which transform differendy under Eq. (51) than i//, and thus would spoil local gauge invariance. One can cure 
this by introducing new vector gauge fields (x) and defining covariant derivatives by 



The gauge variation of this from the variation of ^ alone (i.e. assuming A" constant for the moment) reads 
5xi,D^ Va(x) = is" {x)taDfi va(x) + i dfis" (x) + C^/ij^ {x)e^{x) ta wix) , 



(53) 



(54) 



where we have used Eq. (52). The new term proportional to the structure constants C^^ results from moving the gauge 
variation of yf in Eq. (53) to the left of the matrix gauge field A" (x)ta and is only present in non-abelian gauge theories 
for which the ta do not all commute. The variation dx^Sm of the matter action 5m can then be obtained from Eq. (48), 
generalized to several e", and with d^e'^ (x) substituted by the corresponding first factor of the second term in Eq. (54), 



-J d'xJ&ix) [a^e«(x)+C^/^(x)£^(x) 



where the gauge currents Eq. (50) now read [compare Eqs. (47) and (51)] 



(55) 



(56) 



in terms of the matter Lagrangian ReaUzing now that \dSra/ dA^{x)\ — \d^ra/ ^(^^1 V^)] ( — ita^f) — Jat we see that 
5Sm = SxifSm + [dSm/ dA" (x)] 5 A" (x) vanishes identically if we adopt the gauge transformation 



5A« (x) = a;,e« (x) + C^«X (x)e^(x) , 

for the gauge field A^(x). 

The standard gauge-invariant term for fermions is then given by the matter Lagrange density 

■^m = -W{fD^+m)\if=-\fr{fd^+m)\i/+Ap^, 

where m is the fermion mass matrix. The second equality shows how the matter Lagrangian s] 
quadratic in the fields, Eq. (26), and the fundamental coupUng of the gauge field to the gauge c 
is real and the gauge current is hermitian, J^a = jjiw the gauge fields A^ are also real. 



(57) 



(58) 



Lectures on Astroparticle Physics 



February 2, 2008 



10 



4.3. Gauge Theory of the Electroweak Interaction 



In the electroweak Standard Model the elementary fermions are arranged into three families or generations which 

V; 



here are labeled with the index Each family consists of a left-chiral doublet of leptons, 



, a left-chiral doublet 



of quarks, 



and the corresponding right-chiral singlets Lg, um, and djR. Here, left- and right-chiral is understood 



as in Eq. (31), and each quark species comes in three colors corresponding to the three-dimensional representations of 
the strong interaction gauge group SU (3) whose index is suppressed here. The three known leptons are the electron, 
muon, and tau with their corresponding neutrinos. The three up-type quarks are called up, charm-, and top-quark, and 
the down-type quarks are the down-, strange-, and bottom-quarks. The fermion masses rise steeply with generation 
from about 1 MeV for the first generation to up to ~ 175 GeV for the top-quark whose direct discovery occurred as 
late as 1995 at Fermilab in the USA. 

Note that no right-handed neutrino appears and thus neutrino mass terms of the form vEVR+h.c. (h.c. denotes 
hermitian conjugate here and in the following) are absent in the Standard Model. Implications of recent experimental 
evidence for neutrino masses for modifications of the Standard Model will not be discussed here. To simplify the 



notation we assemble all fields into lepton and quark doublets, /, 
components. We will also use the Pauli matrices 



1 ,and^,: 



including the right-handed 



(tb,T) = ('rb,Ti,T2,T3) 

'I OW 1/2 

lj'U/2 



(59) 



-i/2 
i/2 



1/2 
-1/2 



The electroweak gauge group is given by 



G = SU{2)lxU{1)y, 



(60) 



where the first factor only acts on the left-handed doublets. Denoting the dimensionless coupling constants correspond- 
ing to these two factors with g and g', we write the four generators in the leptonic and quark sector as 



1 + 75 

U^tq = {tut2,t3)^g—r—T: 



tYl = g 



tYq = g' 



1 + 75 To 1 - 75 
1 + 75^0 1-75 /To 



(61) 



These correspond to the generators from the previous section, and we denote the corresponding gauge fields by 
and B^. It is easy to see that the electric charge operator is then given by the combination 



I- 



- f3 jtY, 



(62) 



where e is the (positive) electric charge unit. 

We are here only interested in the part of the Lagrangian involving matter fields. This is then given by Eq. (58) 
where i/a now represents all lepton and quark multiplets /; and^j. Using Eq. (53), where, from comparing Eq. (61) with 
Eq. (52), Cp^ = gEapy for SU (2)/,, and zero for U (1), we can write the matter part of the electroweak Lagrangian as 



- £ liY^ {dn - i\n ■ t - iBntYi) U 

i=l 
3 

- qi'f {dn - iA^ ■ t - iB^tYq) qi ■ 

i=l 



(63) 
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It will be more convenient to use charge eigenstates as basis of the electroweak gauge bosons and to identify the 
photon as carrier of the electromagnetic interactions. There is then one other neutral gauge boson Z^u and two gauge 
bosons of charge ±e. They are defined by 



cos 0ewZ^ - sin 9ewA^ 
= sin 0ewZ;u + cos Q^v/A^i , 



A, = 



Af, = 



(64) 



(65) 



where the electroweak angle 9ew is defined by 

g = 



sin^ew 

The interaction terms in Eq. (63) can then be written as 

g 1+75 



cos dm 



V2 2 

8 



{W+x++W^T-) 



(66) 



(67) 



2 cos 



'75 -2 

; — T3-^sm 0OW 



where T^*" = Ti ± iT2 are the weak isospin raising and lowering operators, respectively. 

Up to this point all fields are massless. Mass terms for gauge bosons and for fermions, whether Dirac or Majorana 
(see below), are inconsistent with gauge invariance. The standard way to introduce them is by spontaneously broken 
gauge symmetries. Without going into any detail here, we just mention that this is done by introducing a scalar Higgs 
field coupUng to gauge boson and fermion biUnears in a gauge-invariant way and making it adopt a vacuum expectation 
value due to a suitably chosen potential. 

Let us now consider processes involving the exchange of a or Z boson with energy-momentum transfer q much 
smaller than the gauge boson mass, |^^| ^ '"wz' ^^'^^ that the boson propagator can be approximated by —iTf^v/fn^z- 
In this case the second order terms in the perturbation series give rise to effective interactions of the form 



— T^^T^ A- — J^J 



2 •'nc-' Unc • 



(68) 



Here, the charged current and neutral current are gauge currents given by comparing Eq. (58) with Eq. (63), and using 
Eq. (67), 



.8 ^ 



1 



•'nc 



8 



3 r 



2 cos ft 



■ew ,=1 



+qi'f 



; — T3 -^Sin 0ew 1 k 



- 75 2 

; — X^-qsm 0ew I qi 



(69) 



Eqs. (68) and (69) provide an effective description of all low energy weak processes. This is an instructive example 
of how a more fundamental renormalizable description of interactions at high energies, in this case electroweak gauge 
theory, can reduce to an effective non-renormalizable description of interactions at low energies which are suppressed 
by a large mass scale, in this case mw or mz- In fact, the latter is identical in form with the historical "V-A" theory 
Eq. (37) which, for example, for the muon decay jii" ^ e~ v<,Vm reads 



Gf 
V2 



e-'f{\ + y5)Ve [v^7m(1 + 75)M ]+h.c., 



(70) 



Lectures on Astroparticle Physics 



February 2, 2008 



12 



where the Fermi constant Gp by comparison with Eqs. (68) and (69) is given by 



o2 

Gf= —^--5- = 1.16637(1) X lO-^GeV-^. (71) 



w 



Radioactive J3— decay processes are described by the terms in Eq. (68) containing e^7^(l + 75)Ve or its hermitian 
conjugate for one of the charged currents J^^c or 7^cc> ™d a quark term for the other current. For example, neutron 
decay, n pe^Ve is due to the contribution uY^{l + Y5)d to 7^, which causes one of the d-quarks in the neutron to 
transform into a u-quark under emission of a boson which in turn decays into e^Vg, {add) {uud)e^Ve. 

Inverting Eq. (64) to = cos 0ew^^ + sin^ewS^ and writing out the mass term of the neutral gauge boson sector 
jM^Z^Z^^ imphes 

niw = mz cos > (72) 

because the mass term of has to be identical to the one for W^. From this it follows immediately that the neutral 
current part of Eqs. (68), (69) involving neutrinos can be written as 

^ [Vifil + r5)Vi] [vYn {8l{1 + 75) +^«(1 - 75)) V] , (73) 



where xj/ stands for quarks and leptons and 



8L = T3-q'sin2 0ew 

gR = -qsin^e^^. (74) 



5. NEUTRINO SCATTERING 

Imagine a neutrino of energy Ey scattering on a parton ; carrying a fraction x of the 4-momentum P of a state X of 
mass M. Denoting the fractional recoil energy of X hy y = E'^/ Ey and the distribution of parton type ( by fi (x, Q), in 
the relativistic limit Ey mx the contribution to the vX cross section turns out to be 

Jd^x 2GlMEyx ( y 

= [ iMEyxy + Ml^, ) [Si,L + 8iA^-y) ] • (75) 

Here, g,- ^ and gi^R are the left- and right-chiral couphngs of parton respectively, given by Eq. (74). Eq. (75) applies to 
both charged and neutral currents, as well as to the case where X represents an elementary particle such as the electron, 
in which case fi{x,Q) = 5 (x — 1 ) . 

As usual, if the four-momentum transfer Q becomes comparable to the electroweak scale, ^ '"wz' weak 
gauge boson propagator effects, represented by the factor 2/ (Q^ + ^wz) ™ ^1- (75), become important. We have 
used that in the hmit \Q^\:$>M^ one has ~ -M^ = {xP + Qf- ~ + 2P ■ Qx with P • 2 ~ -ME'^ = -ME^y evaluated 
in the laboratory frame, i.e. the rest frame of X before the interaction. ~ IMEyxy is also called the virtuality because 
it is a measure for how far the exchanged gauge boson is form the mass shell = -M2. 

We will not derive Eq. (75) in detail, but it is easy to understand its structure: First, the overall normalization is 
analogous to Eq. (8), using the fact that for Ey^M the CM momentum p^, ~ MEy /2. Second, if the helicities of the 
parton and the neutrino are equal, the total spin is zero and the scattering is spherically symmetric in the CM frame. 
In contrast, if the parton is right-handed, the total spin is 1 which introduces an angular dependence: After a rotation 
by the scattering angle 0« in the CM frame the particle helicities are unchanged for the outgoing final state particle 
and one has to project back onto the original hehcities in order to conserve spin. If a left-handed particle originally 
propagated along the positive z-axis, its left-handed component after scattering by 0* in the x — z plane is 

(76) 



2 v p J \ ^ J 2 V i+cosa 

giving a projection [( 1 + cos 0* )/2]^. Now, Lorentz transformation from the CM frame to the lab frame gives Ey /Ey = 
(1 + vxcos0*)/2 ~ (1 -|-cos0*)/2 in the relativistic limit and thus the projection factor equals {Ey/Ey)^ = (1 — 
E'^/Eyf- — (1 —y)^, as in Eq. (75) for the right-handed parton contribution. Integrated over < y < 1 this gives 1/3, 
corresponding to the fact that only one of the three projections of the 7=1 state contributes. 
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\ [GeV] 

FIGURE 1. FromRef. [5]. Cross sections for AT interactions at high energies, according to the CTEQ4-DIS parton distributions: 
dashed line, a{veN Ve + anything); thin line, a{V(N — » £~ + anything); thick Une, total (charged-current plus neuttal-current) 
cross section. 



5.1. Neutrino-Nucleon Scattering and Applications 

We now briefly consider neutrino-nucleon interaction. From Eq. (75) it is obvious that at ultra-high energies 
2EyM ^ Myi/z^ the dominant contribution comes from partons with 

Since, very roughly, xfi{x,Q) oc x^"-^ for x <C 1, it follows that the neutrino-nucleon cross section grows roughly 
oc Ey^. This is confirmed by a more detailed evaluation of Eq. (75) shown in Fig. 1. 

Let us use this to do a very rough estimate of event rates expected for extraterrestrial ultra-high energy (UHE) 
neutrinos in neutrino telescopes. Such neutrinos are usually produced via pion production by accelerated UHE protons 
interacting within their source or with the cosmic microwave background (CMB) during propagation to Earth. The 
threshold for the reaction Ny Nk, for a head-on colUsion of a nucleon N of energy E with a photon of energy e is 

given by the condition s={E + e)'^ — [(E^ — rn^yi'^ — e] > (otat + w^t)^, or 

'^>«^-3.4xlO"(^)-'eV. ,78, 

where e ~ 10^^ eV represents the energy of a typical CMB photon. At these energies, the secondary neutrino flux 
should therefore be very roughly comparable with the primary UHE cosmic ray flux, within large margins. Fig. 2 
shows a scenario where neutrinos are produced by the primary cosmic ray interactions with the CMB. Using that the 
neutrino-nucleon cross section from Fig. 1 roughly scales as Oyn °< Ey^^^ for 10^^ eV <Ev< 10^^ eV, and assuming 
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E[eV] 

FIGURE 2. A scenario from Ref. [6, 7] where neutrinos are produced as secondaries of UHE cosmic rays (data with error bars 
above lO'^ eV) by interactions with the cosmic microwave background (CMB). Theoretical fluxes of protons, 7— rays, and neutrinos 
(per flavor) are as indicated. Also shown are the atmospheric neutrino flux [8], as well as existing upper limits on the diffuse neutrino 
fluxes from MACRO [9], AMANDA II [10], BAIKAL [11], AGASA [12], the Fly's Eye [13] and RICE [14] experiments, and the 
limits obtained with the Goldstone radio telescope (GLUE) [15] and the FORTE satellite [16], as indicated. The cosmic ray data 
are from the AGASA [17] and HiRes [18] experiments To the left the newest estimate of the diffuse GeV y— ray background from 
EGRET data is shown [19]. Since most of the electromagnetic energy ends up in this energy range and since pion production 
produces comparable amounts of y— rays and neutrinos, predicted neutrino fluxes cannot considerably overshoot the horizontal 
line on the level of the EGRET estimate marked y— ray bound. If the sources are transparent to produced cosmic rays, the more 
restrictive but less general "Waxman-Bahcall bound" [20], marked "WB bound", results. This will be discussed in more detail in 
Sect. 9.2. 



water or ice as detector medium, we obtain the rate 

Tv ~ avNiEv)27tEvjiEv)nNVef[ 

{lO^'^ewJ lvl02eVcm-2sr-is-V UmV ' ^ ' 

where ~ 6 x lO^^cm^^ is the nucleon density in water/ice, Veff the effective detection volume, and j{Ev) is the 
differential neutrino flux in units of cm^^e ' sr^ ^ s^ ' . 

Eq. (79) indicates that at Ey > lO'^eV, effective volumes > 100 km^ are necessary. Although impractical for con- 
ventional neutrino telescopes, big air shower arrays such as the Pierre Auger experiment can achieve this. In contrast, 
if there are sources such as active galactic nuclei emitting at Ey ^ lO^^eV at a level E^ j{Ev) ^ lO^eVcm^^sr^'s^^ 
km-scale neutrino telescopes should detect something. Such fluxes are consistent with general considerations, see 
Fig. 2. 

Finally, Fig. 3 shows more detailed neutrino flux sensitivities expected from future experiments. 
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FIGURE 3. From Ref. [6]. Similar to Fig. 2, but showing expected sensitivities of the currently being constructed Pierre Auger 
project to tau-neutrinos [21], the planned projects Telescope Array (TA) [22], the fluorescence/Cerenkov detector NUTEL [23], the 
space based EUSO [24], the water-based Baikal [11] and ANTARES [25] (the NESTOR sensitivity for 1 tower would be similar 
to AMANDA-II and for 7 towers similar to ANTARES [26]), the ice-based AMANDA-II [10] and ICECUBE [27] (similar to the 
intended Mediterranean km-' project [28]), and the radio detectors RICE [29] and ANITA [30], as indicated. All sensitivities except 
for ANITA and RICE refer to one year running time. For comparison, the 7— ray bound derived from the EGRET GeV 7— ray 
flux [19] is also shown. 



6. DIRAC AND MAJORANA NEUTRINOS 

Up to now we have assumed that neutrinos and anti-neutrinos are separate entities. This is true if lepton number is 
conserved, see Tab. 1, and corresponds to pure "Dirac neutrinos". However, lepton number may be violated in the 
neutrino sector and neutrinos may be indistinguishable from anti-neutrinos. In order to elucidate this, let us first study 
some symmetries of the Dirac equation (27). From Eqs. (23), (25) one can easily show that 

1^= 72 7m 72- (80) 

Complex conjugating the Dirac equation (27) and multiplying it with £,*Y2 from the left, it then follows that it is 
invariant under the "charge conjugation transformation" 

Cxi/{x)C-' = w' = ^*r2W*, (81) 

where £, is an arbitrary complex number with |^| = 1. This transformation exchanges particles and anti-particles 
and satisfies = y/. Note that 72 appears because according to Eq. (80) it is the only real Dirac matrix in this 

convention. 

A Majorana neutrino satisfies the reality condition 

(/)(x) = 72(/)*(x), or r = ^*<P- (82) 

These spinors are of the form 
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where ;|f is a two-spinor. This implies that for any Dirac spinor \j/ one can construct a Majorana spinor by 



= (84) 

Note that this spinor is not an eigenstate of lepton number because under a phase transformation ^ — > ^e'" one has 
^e~'". Defining 

wi,R ^ {r\R = ^-^r = (^-^wY , (85) 



2 

one can define left- and right-handed Majorana fields 

<l>± = Wlm + iWLfiY = WLji + Wr,l ■ (86) 

Note that both these fields now contain both left and right-handed fields. What before experimentally was called 
neutrino and anti-neutrino now is called left- and right-handed neutrino, respectively. We can now introduce Majorana 
mass terms of the form 

-^M = ~{mi^<i)++m]^<i)-) (87) 

= -\mL [WWr + V^V^l) - \mR {WVl + Wlw) , 
where m/, and iur are real. Together with the Dirac mass term this can be written as 

where we have used (v/T'Ki)^ = W2WI' around Tab. 2 and (using 72 = 72 and the fact that y/ anti-commutes) 



WlWr 



for any Dirac spinors 2 and \\f. Note that under \\f — > \\fe^" Dirac terms are invariant, whereas Majorana terms 
pick up the phase e*^'", according to lepton number conservation and non-conservation, respectively. Furthermore, 
we see that for niR ^ mo, niL — 0, the two mass eigenvalues in Eq. (88) are ~ niR and m^/niR. The latter are very 
small and thus may explain the sub-eV masses involved in left-chiral neutrino oscillations. This is called the see-saw 
mechanism which would imply that the mass eigenstates are Majorana in nature. The existence of one heavy right- 
handed Majorana neutrino per lepton generation is motivated by Grand Unification extensions of the electroweak 
gauge group to 50(10) which has 16-dimensional representations that could fit 15 Standard Model lepton and quark 
states plus one new state, see, e.g., Ref. [31]. 

Finally, in the exactly massless case, Dirac and Majorana particles are exactly equivalent, since the two fields 
Eq. (86) completely decouple, see Eq. (88). We also mention that in supersymmetric extensions of the Standard Model, 
the fermionic super-partners of the gauge bosons are Majorana fermions. As a consequence, they can self-annihilate 
which plays an important role to their being candidates for cold dark matter. 

For n > 1 neutrino flavors, mass eigenstates |v,) of mass m,- and interaction eigenstates |Va) in general are not 
identical, but related by a unitary nxn matrix U : 

\Va) =Y,Uai\Vi) , (89) 

i 

where for anti-neutrinos U has to be replaced by U*. Such a matrix in general has real parameters. Subtracting 
2n — 1 relative phases of the n neutrinos in the two bases, one ends up with (n — 1)^ physically independent real 
parameters. Of these, n(n— l)/2 are mixing angles, and the remaining (n — l)(n — 2)/2 are "CP- violating phases". In 
order to have CP- violation in the Dirac neutrino sector thus requires n > 3. Once the relative phases of the different 
flavors have been fixed, for non- vanishing Majorana masses there will in general be n — 1 Majorana phases that can not 
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be projected out hy y/^ y/e'" in Eq. (87). Thus, the number of independent real parameters is larger, namely n{n — 1), 
in this case. Note that the corresponding Cabibbo Kobayashi Maskawa ( CKM) matrix in the quark sector is pure Dirac 
because Majorana terms would violate electric charge conservation in the quark sector. 

If at time r = a flavor eigenstate \Va) = T.iUai\Vi) is produced in an interaction, in vacuum the time development 
will thus be 

|V(0> = Y^Uate-'^'' \Vi) = '£UaiU*p,e-'^-' \v^) . (90) 

Since masses and energies of anti-particles are equal according to the CPT theorem, from this we obtain the foUowing 
transition probabiUties 



£t/„,f/^*,.exp(-/£,f) 

i 



(91) 



From this foUows immediately 

PiVa ^ Vp) = P{Vp ^ Va) , (92) 
which is due to the CPT theorem. Furthermore, if the mixing matrix satisfies a reality condition of the form 

Uai = Kirii, (93) 

with r]i phases, corresponding to Cf -conservation, one also has 

P{Va ^ Vp) = P{Va ^ Vp) . (94) 

We mention two other important differences between Dirac and Majorana neutrinos: 

• Neutrino-less double beta-decay is only possible in the presence of Majorana masses because the final state e~e~ 
violates lepton number. In this case the rate is proportional to the square of 



rUee = 



(95) 



see Eq. (88), where only one of the Majorana phases a, can be projected out. Apart from these phases, this 
equation results from Eq. (90) for a = J3 = e. There is even evidence claimed for this kind of decay, and thus 
for an electron neutrino Majorana mass around 0.4 eV, see Ref. [32]. The issue is expected to be settled by next 
generation experiments such as CUORE [33]. 

In contrast, in j3 —decay with neutrinos, the electron spectra are influenced by the individual eigenstates of real 
mass Mi, and not by any phases. The current best experimental upper limit is given by the Mainz experiment 
based on tritium J3— decay — >^ Hee~ v,, [34], 



(96) 



at 95% confidence level (CL). The KATRIN experiment [35] aims at a sensitivity down to 0.2 eV within the next 

few years. 

Majorana neutrinos cannot have magnetic dipole moments between equal neutrino flavors, as seen from the 
following identity using Eqs. (82), (80), (34), and the reality of the spinors in Tab. 2: 

^Wl[7^l,7vW2 - -wI72f[7^l,7v]r2¥2 = -Wi (r^fe, 7v])* Wi = 



7vm 



-iV2[7^i,7v]Wl: 



where in the second-last identity we have used that transposition changes the order of the fermionic fields, thus 
picking up a minus sign. As a consequence, only transition magnetic moments between different flavors are 
possible for Majorana neutrinos. 



Lectures on Astroparticle Physics 



February 2, 2008 



18 



7. NEUTRINO OSCILLATIONS 



Let us now restrict to two-neutrino oscillations, n = 2, between | Ve ) and \ vn), say, and write 

jj _ ( '-OS ^0 sin 00 
sin 00 cos 00 



(97) 



for the mixing matrix in Eq. (89) which is characterized by one real vacuum mixing angle 0o. Since id | v,) /dt = Ei | V;) 
for J = 1 , 2 in the mass basis, and since = (m? + p^) p+mf / {2p) c±E+m]/ {2E) in the relativistic limit p ^ nij, 
using the trigonometric identities cos^Qq — sin^Qo = cos20o. 2cos0osin0o = sin20o, it follows from Eqs. (89), (97) 
that 



dt 



E+- 



m\+m^\ Am^ / cos20o 



AE 



+ 



-sin 2 00 
AE V-sin20o -cos20o 



(98) 



where we consider a given momentum mode p, and Am^ = m\ — m\. From now on we will consider the relativistic 
limit with p c:iE. The first term in Eq. (98) is a common phase factor and can be ignored. 

The integrated version of this is Eq. (90). Then applying Eq. (91), one can show that this has the solution 



P{Va ^ Vjs) = i sin2 20o - cosAw^ for a ^ J3 , 



for oscillations over a length L. The oscillation length in vacuum is thus 



Lo = An 



jAm^ 



:2.48 



VMeV 



|Am2 
eV^ 



m. 



(99) 



(100) 



Neutrino oscillations are modified by forward scattering amplitudes in matter. Since neutral currents are by definition 
flavor-neutral, they only contribute to the common phase factor which in the following will be ignored. The charged 
current interaction is diagonal in flavor space and, according to Eqs. (68), (69), and (71), the low-energy limit of its 
forward scattering part for Ve(p) has the form 



2^/2Gp£ L(p')t70/i±^v,(p)l L(p)t/r^i±^.(p' 



(101) 



where we have used y/ = (//^/y". We need to express this in terms of Dirac bilinears of the form Tab. 2 for electrons 
and neutrinos separately. In order to do that we use the fact that every 4x4 matrix can be expanded according to 



' = 1 



tr 



tri 



(102) 



where are the 16 matrices appearing in Tab. 2 which satisfy Xv{0'iffj) = for i ^ j. Using this one can show that 



'ft 



1 + 75 



1+75 



1+75 

2 



r/3 



a/3 
1 + 75 



aS 



yS 



1 + 75 



rP 



a8 



(103) 



Applying this to Eq. (101) and noting that the last term in Eq. (103) does not contribute in the rest frame of the electron 
plasma, one obtains 

-2\/2Gf 



v.(p)^-^v,(p) 



e(p')^^e(p') 



(104) 



where we have picked up an extra minus from the anti-commutation of fermionic fields. The electron and neu- 
trino fields have the form Eq. (17). The matter-dependent part of the sum in Eq. (104) thus takes the form 
£p(aJ^(p)ae,L(p) — fog2,(P)^e,L(p)- When tracing out the charged lepton density matrix, this reduces to the density 
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of left-chiral electrons minus the density of left-chiral positrons. Thus, for an unpolarized plasma, the contribution to 

the \Ve) self-energy finally is —^/iGpNe, where Ne is the electron-number density, i.e. the electron minus the positron 
density, and analogously for the other active flavors. The non-trivial part of Eq. (98) is thus modified to 



Am cos20Q /o/^ Am sin20Q 

( I ^1 4£ V^<-rF^Ve I ( 1 (105) 



4E 




It is illustrative to write this in terms of the hermitian density matrix pp{t) = |Vp(f)) (Vp(f) |. If we expand this into 
an occupation number np and polarization Pp, Pp = j (np -|-Pp- cj), one can easily show that Eq. (105) is equivalent to 



Pp=BpXPp, (106) 



where the precession vector Bp = B(p) 



Blip) = ^sm20o, 

Hp) = ^cosIOo-ViGfN, (107) 

where N = Ng — N^^^ for Vg — V^^^ mixing and = A^e,;u,T — N„ for active-sterile mixing, with N„ a combination of 
nucleon densities. Thus, neutrino oscillations are mathematically equivalent to the precession of a magnetic moment 
in a variable external magnetic field. 

Eq. 107) shows immediately that at the resonance density 

^ _ Am^ cos 200 ^ 

2V2GfE ^^^^^ 

the two diagonal entries in Eq. (105) become equal or, equivalently, B^{p) = 0. More generally, Eqs. (105), (106), 
(107) are diagonalized by a mixing matrix Eq. (97) where 0o is replaced by the mixing angle in matter 9e given by 

Maximum mixing, 9E = n/4, thus occurs at the so-called Michaev-Smirnow-Wolfenstein (MSW) resonance at N = Nr. 
For N Nr one thus has vacuum mixing, Oe — Oq, whereas for N ^ Nr one has Oe «C 9o- Eq. (106) now shows 
that propagation from N ^ Nr to N ^ Nr can lead to an efficient transition from one flavor to another, as long as 
I Bp I < I Bp p. Such transitions are called adiabatic. Note that since masses of anti-particles and particles are equal, 
whereas the lepton numbers Ng^^^x change sign under charge conjugation, resonances in matter occur either for 
neutrinos or for anti-neutrinos. 
In subsequent sections we will apply this to various observational evidence for neutrino osciUations. 



8. SELECTED APPLICATIONS IN ASTROPHYSICS AND COSMOLOGY 
8.1. Stellar Burning and Solar Neutrino Oscillations 

Weak interactions are crucial in cosmology and steUar physics. In main sequence stars the first stage of hydrogen 
fusion into helium is the weak interaction 

p + p^^U + e+ + Ve. (110) 
The subsequent reactions ^Yl+p — >^He-|-7 and •'He-l-^He— >^He-|-p -|- p lead to the net reaction 

4p^'^Ue + 2e+ + 2Ve + 26J3MeW. (Ill) 

For many more details see Ref. [36]. 
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FIGURE 4. The predicted solar neutrino energy spectrum from Ref. [37]. For continuum sources, the neutrino fluxes are given in 
number per cm^^scc^'MeV^' at the Earth's surface. For line sources, the units are number per cm~^sec~^. The total theoretical 
uncertainties are shown. The CNO neutrino fluxes are not shown. 



When normalized to the solar energy flux arriving at Earth, one can calculate the expected neutrino fluxes within 
the so-called Standard Solar Model. The resulting fluxes in this "pp" channel as well as in other reaction channels are 
shown in Fig. 4 

However, less than half of the expected solar electron-neutrino flux at a few MeV has been observed. On the other 
hand, neutral current experiments with the Sudbury Neutrino Observatory (SNO) [38] have shown that the sum of the 
electron, muon- and tau neutrino flux coincides with the expected electron neutrino flux. This can be explained by 
an MSW transition of Vg into v^, and within the Sun with a Am^^j^ ~ 10~^eV^. Note that this corresponds to a 
vacuum oscillation length Eq. (100) of a few hundred kilometers. Recently this has been confirmed independently by 
the KamLAND experiment [39] which measured the disappearance of the Vg neutrinos produced by nuclear reactors 
a few hundred kilometers from the detector. The best fit parameters for the parameters of mixing of two neutrinos in 
vacuum from aU solar and reactor data are [40] 



Aml,^ ^ S.2tr,{t'oi) X lO-^eV^ ; tan^ ^ 039t'oXil) , 



(112) 



where 1 a and 3ct errors are given. The relevant contour plots are shown in Fig. 5. It is interesting to note that maximal 
mixing is strongly excluded, tan 0soiar < 1-0 at 5.8(7. 

Finally, since solar and reactor neutrino experiments deal with Vg and Vg, respectively, comparison of the two 
corresponding osciUations parameters allows to set limits on CPT-violation in the neutrino sector. 
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FIGURE 5. From Ref. [40]. Allowed oscillation parameters: Solar vs KamLAND. The two left panels show the 90%, 95%, 99%, 
and 3cT allowed regions for oscillation parameters that are obtained by a global fit of all the available solar data. The two right 
panels show the 90%, 95%, 99%, and 3cT allowed regions for oscillation parameters that are obtained by a global fit of all the 
reactor data from KamLAND and CHOOZ. The two upper (lower) panels correspond to the analysis of all data available before 
(after) the Neutrino 2004 conference, June 14-19, 2004 (Paris). The new KamLAND data are sufficiently precise that matter effects 
discernibly break the degeneracy between the two mirror vacuum solutions in the lower right panel. 
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FIGURE 6. From Ref. [41]. Allowed (sin^ 0atmos, AOT^t^no^) regions at 90%, 95%, 99%, and 3a CL for two degrees of freedom. 
The regions delimited by the lines correspond to atmospheric data only, while for the colored regions also K2K data are added. 
The best fit point of atmospheric (atmospheric + K2K) data is marked by a triangle (star). Also shown is the Ax^ as a function of 
sin^ Oatmos and Am^i^^^^, minimized with respect to the undisplayed parameter. 



8.2. Atmospheric Neutrinos 

Cosmic rays interact in the atmosphere and produce, among other particles, pions and kaons whose decay products 
contain neutrinos. The observed ratio of upcoming to down-going atmospheric muon neutrinos is about 0.5 in the GeV 
range. Since upcoming neutrinos travel several thousand kilometers, this can be interpreted as vacuum oscillations 
between muon and tau-neutrinos. The best fit parameters are 

Am2j^„,~ 2.1 X 10-3 eV2; sin^ e^tmos ^ 0.5 , (113) 

consistent with maximal mixing. The corresponding contours resulting from atmospheric and long baseline neutrino 
oscillation data from K2K [42] are shown in Fig. 6. Recently, the L/E dependence, see Eq. (99), characteristic for 
neutrino oscillations has been confirmed by the Superkamiokande experiment [43], thereby strongly constraining 
alternative explanations of muon neutrino disappearance such as neutrino decay [44]. In addition, oscillations into 
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sterile neutrinos are strongly disfavored over oscillations into tau-neutrinos via the following discriminating effects: 
Neutral currents would be non-diagonal for oscillations into sterile states, thus modifying oscillation ampUtudes and 
total scattering rates, and charged current interactions of tau-neutrinos imply T appearance. 

In 3-neutrino oscillation schemes ©solar and ©atmos are usually identified with du and 623, respectively (we here 
assume the "normal mass hierarchy" mi , 1112 <C W3 which seems the most natural for neutrino mass modeling in Grand 
Unification scenarios [45]). According to the discussion around Eq. (89), there is one more mixing angle 613 and 
at least one Dirac CP- violating phase called 5. Note that solar and atmospheric neutrinos only decouple exactly for 
013 = in which case CP would also be conserved. Whereas there are at most weak indications for leptonic CP- 
violations yet [46], the third mixing angle is constrained at 3c7 CL by Ref. [41] 

sin^ 013 < 0.061. (114) 

We finally stress that neutrino oscillations are sensitive only to differences of squared masses, not to absolute 
mass scales. To probe the latter requires laboratory experiments discussed earlier such as j3— decay, the study of 
cosmological effects such as the influence of neutrino mass on the power spectrum, see Sect. 8.4, or measuring time 
delays of astrophysical neutrino bursts from 7— ray bursts and supemovae relative to the speed of light. 



8.3. Big Bang Nucleosynthesis (BBN) 

For more detailed introductions to the following three topics we refer the reader to standard text books [47, 3 1]. 

The early universe consisted of a mixture of protons, neutrons, electrons, positrons, photons and neutrinos. Their 
relative abundances were determined by thermodynamic equilibrium until the weak inlcraclions "froze out" once the 
temperature of the expanding universe dropped below 7/ ~ 1 MeV where their rates became smaller than the expansion 
rate. For example, according to Eqs. (7) and (8) the interaction rates of nucleons VeP ne~^ and e~p nVg are 

r^naocGlT^ (115) 

at temperatures 100 GeV >T> IMeV where the neutron-proton mass difference m„—mp = 1.293 MeV and the 
electron mass are negligible and the and electron neutrino densities n^T^. This becomes indeed compjirable to 
the expansion rate 

H'^p^l^/ni'e^gV^T^/mv, (116) 

where p is the total energy density and the number of relativistic degrees of freedom, once T approaches Tf ~ 
1 MeV. The equiUbrium neutron to proton ratio at that temperature is given by thermodynamics as 

— =exp[-(OT„-OTp)/7y] (117) 

Hp 

At that time, the free neutrons were quickly bound into helium which could not be broken up any more by the cooUng 
thermal radiation. The helium abundance was thus determined by the freeze out of electroweak interactions. Since 

equating Eq. (115) with Eq. (116) yields Tj oc we also see that the hehum abundance should increase with g*. 
Since the number of stable neutrino species with mass below ~ 1 MeV contributes to g*, this number is constrained 
by the observed helium abundance. More generally, in the absence of a significant asymmetry between neutrinos and 
anti-neutrinos, elemental abundances depend only on the effective number of relativistic neutrinos and the baryon 
to photon ratio 

r/io^lQio^^. (118) 

Hy 

Predictions for standard big bang nucleosynthesis (SBBN) with A^v = 3, the number of active neutrinos consistent with 
the Z boson width, are shown in Fig. 7. A detailed comparison of measured and predicted abundances shown in Fig. 7 
with 7710 and A^v free parameters yields the following: The universal density of baryons Tjio inferred from SBBN and 
the measured deuterium abundance, 7710 (SBBN) = 6.10^q in excellent agreement with the baryon density derived 
largely from CMB data [48], T7io(CMB) = 6.14 ±0.25. However, there is a ~ 2a tension between the ^He abundance 
predicted by SBBN with this concordance Tjio and the observed one. This tension can be mitigated if A^v is allowed to 
be smaller than the canonical A^v — 3. If both the baryon density tjiq and A^v, or equivalently, the expansion rate, are 
aUowed to be free parameters, BBN (D, ^He, and ^^He) and flie CMB (WMAP) agree at 95% CL for 5.5 < tJiq < 6.8 
(0.020 < Q.^h^ < 0.025 for the baryon density in terms of the critical density) and 1.65 < Ny < 3.03 [49]. Are these 
hints for new physics ? 
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FIGURE 7. From Ref. [49] . The SBBN-predicted primordial abundances of D, ^He, and ^Li by number with respect to hydrogen, 
and the ^He mass fraction Yp, as a function of the nucleon (baryon) abundance parameter rjiQ. The bands reflect the theoretical 
uncertainties (±1ct) in the BBN predictions. 



8.4. Neutrino Hot Dark Matter 

Finally, massive neutrinos in the eV range contribute to the density of non-relativistic matter in today's universe, 

n,/z2 = ^!!!^ (119) 

92.5eV ^ ' 

in terms of the critical (closure) density for niy ^ 10"^ eV, today's temperature. Eq. (119) results from the fact that 
neutrinos have been relativistic at decoupling at T ^ 1 Me V, thus constituting hot dark matter, and their number density 
is simply determined by the redshifted number density at freeze-out, in analogy to Sect. 8.3. 

Since neutrinos are freely streaming on scales of many Mpc, the matter power spectrum is reduced by a relative 
amount AP,„/P,„ = — Sily /^m, where i2„, is the total matter density. A combination of data on the large scale structure 
and the CMB then leads to the limit [50] 

rmv<1.0eV. (120) 
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There was even a claim for a positive detection with [5 1 J 

Y^myC^O.SeeV, (121) 

but newest analyses suggest upper bounds even shghtly below this [52]. 

It is intriguing that direct experimental bounds Eq. (96) and cosmological bounds Eq. (120) have reached compara- 
ble sensitivities. In addition, both a combination of future CMB data from the Planck satellite with large scale structure 
surveys [50] and next generation laboratory experiments such as KATRIN will probe the 0.1 eV regime. 

Assuming three active neutrino oscillations with the parameters discussed in Sects. 8.1 and 8.2 has an interesting 
cosmological consequence: Flavor equilibrium is reached before the BBN epoch and the asymmetry parameter 
i^v = Mv/r, where Hv is the common neutrino chemical potential, is constrained by [53] 

|'^v|<0.07. (122) 

As a consequence, neutrino degeneracy is unobservable in the large scale structure and the CMB. 



8.5. Leptogenesis and Baryogenesis 

Neutrino masses may also play a key role in explaining the fact that we Uve in a universe dominated by matter rather 
than anti-matter The heavy right-handed Majorana neutrinos involved in the seesaw mechanism discussed in Sect. 6 
could have been produced in the early Universe and their out-of-equilibrium decays could give rise to a non- vanishing 
net lepton number L. Non-perturbative quantum effects related to the non-abelian character of the electroweak 
interactions can translate this into a net baryon number B while conserving B — L. The amount of baryon number 
riB created in this scenario is related to the low-energy leptonic Cf- violation phase 5 [54]. Its compatibility with the 
observed value for the baryon per photon number ng/n^ c± 6 x 10~ ^'^ impUes a lower bound mr > few lO^** GeV. Via the 
see-saw relation my — /mr for the light neutrino mass, this corresponds to an optimal range 10-3 eV< my <0.1eV, 
in remarkable agreement with the observed atmospheric and solar neutrino mass scales [54]. In general baryogenesis 
requires violation of baryon number B, charge conjugation C, combined charge and parity conjugation CP, and a 
departure from thermal equilibrium, usually caused by the expansion of the Universe. These conditions are known as 
the Sakharov conditions [55]. For more details see Refs. [31] and [47]. 



9. ULTRA-HIGH ENERGY COSMIC RADIATION 

In the final part we discuss some current theoretical issues around ultra-high energy cosmic rays, 7— rays and neutrinos. 
We will see how some of the topics discussed in the previous two parts play an important role in this subject. 



9.1. Introduction 

High energy cosmic ray (CR) particles are shielded by Earth's atmosphere and reveal their existence on the ground 
only by indirect effects such as ionization and showers of secondary charged particles covering areas up to many 
km^ for the highest energy particles. In fact, in 1912 Victor Hess discovered CRs by measuring ionization from a 
balloon [56], and in 1938 Pierre Auger proved the existence of extensive air showers (EAS) caused by primary particles 
with energies above 10^^ eV by simultaneously observing the arrival of secondary particles in Geiger counters many 
meters apart [57]. 

After almost 90 years of research, the origin of cosmic rays is still an open question, with a degree of uncertainty 
increasing with energy [58]: Only below 100 MeV kinetic energy, where the solar wind shields protons coming from 
outside the solar system, the sun must give rise to the observed proton flux. Above that energy the CR spectrum exhibits 
little structure and is approximated by broken power laws °^ E^^: At the energy £ ~ 4 x lO'^ eV called the "knee", the 
flux of particles per area, time, solid angle, and energy steepens from a power law index 7~ 2.7 to one of index ~ 3.0. 
The bulk of the CRs up to at least that energy is beUeved to originate within the Milky Way Galaxy, typically by shock 
acceleration in supernova remnants. The spectrum continues with a further steepening to 7 ~ 3.3 at £ ~ 4 x 10^^ eV, 
sometimes called the "second knee". There are experimental indications that the chemical composition changes from 
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FIGURE 8. From Ref. [62]. The cosmic ray all particle spectrum. 



light, mostly protons, at the knee to domination by iron and even heavier nuclei at the second knee [59]. This is in fact 
expected in any scenario where acceleration and propagation is due to magnetic fields whose effects only depend on 
rigidity, the ratio of charge to rest mass, Z/A. This is true as long as energy losses and interaction effects, which in 
general depend on Z and A separately, are small, as is the case in the Galaxy, in contrast to extra-galactic cosmic ray 
propagation at ultra-high energy. Above the so called "ankle" or "dip" at E 2:^5 x lO'^ eV, the spectrum flattens again to 
a power law of index y ~ 2.8. This latter feature is often interpreted as a cross over from a steeper Galactic component, 
which above the ankle cannot be confined by the Galactic magnetic field, to a harder component of extragalactic 
origin. The dip at £ ~ 5 x lO'^eV could also be partially due to pair production by extra-galactic protons, especially 
if the extra-galactic component already starts to dominate below the ankle, for example, around the second-knee [60]. 
This latter possibility appears, however, less likely in light of a rather heavy composition up to the ankle suggested 
by several experiments [59]. In any case, an eventual cross over to an extra-galactic component is also in fine with 
experimental indications for a chemical composition becoming again lighter above the ankle, although a significant 
heavy component is not excluded and the inferred chemical composition above ^ 10'**eV is sensitive to the model of 
air shower interactions and consequently uncertain presently [61]. In the following we will restrict our discussion on 
ultra-high energy cosmic rays (UHECRs) above the ankle where the spectrum seems to continue up to several hundred 
EeV (1 EeV= 10^^ eV) [17, 18], corresponding to about 50 Joules. The all-particle spectrum is shown in Fig. 8. 

We note that until the 1950s the energies achieved with experiments at accelerators were lagging behind observed 
CR energies which explains why many elementary particles such as the positron, the muon, and the pion were first 
discovered in CRs [63]. Today, where the center of mass (CM) energies observed in collisions with atmospheric nuclei 
reach up to a PeV, we have again a similar situation. In addition, CR interactions in the atmosphere predominantly 
occur in the extreme forward direction which allows to probe non-perturbative effects of the strong interaction. This 
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is complementary to collider experiments where the detectors can only see interactions with significant transverse 
momentum transfer. 

Although statistically meaningful information about the UHECR energy spectrum and arrival direction distribution 
has been accumulated, no conclusive picture for the nature and distribution of the sources emerges naturally from 
the data. There is on the one hand the approximate isotropic arrival direction distribution [64] which indicates that 
we are observing a large number of weak or distant sources. On the other hand, there are also indications which 
point more towards a small number of local and therefore bright sources, especially at the highest energies: First, 
the AGASA ground array claims statistically significant multi-plets of events from the same directions within a few 
degrees [65, 64], although this is controversial [66] and has not been seen so far by the fluorescence experiment 
HiRes [67]. The spectrum of this clustered component is °^ E~^-^ and thus much harder than the total spectrum [65]. 
Second, nucleons above ~ 70EeV suffer heavy energy losses due to photo-pion production on the cosmic microwave 
background — the Greisen-Zatsepin-Kuzmin (GZK) effect [68] already mentioned in Sect. 5.1 — which limits 
the distance to possible sources to less than ~ 100 Mpc [69]. For a uniform source distribution this would predict 
a "GZK cutoff", a drop in the spectrum. However, the existence of this "cutoff" is not estabUshed yet from the 
observations [70] and may even depend on the part of the sky one is looking at: The "cutoff' could be mitigated in the 
northern hemisphere where more nearby accelerators related to the local supercluster can be expected. Apart from the 
SUGAR array which was active from 1968 until 1979 in Austraha, all UHECR detectors completed up to the present 
were situated in the northern hemisphere. Nevertheless the situation is unclear even there: Whereas a cut-off seems 
consistent with the few events above 10^" eV recorded by the fluorescence detector HiRes [18], it is not compatible 
with the 11 events above 10^*^ eV measured by the AGASA ground array [17]. It can be remarked, however, that 
analysis of data based on a single fluorescence telescope, the so-called monocular mode in which most of the HiRes 
data were obtained, is complicated due to atmospheric conditions varying from event to event [71]. The solution of 
this problem may have to await more analysis and, in particular, the completion of the Pierre Auger project [72] which 
will combine the two complementary detection techniques adopted by the aforementioned experiments and whose 
southern site is currently in construction in Argentina. 

This currently unclear experimental situation could easily be solved if it would be possible to follow the UHECR 
trajectories backwards to their sources. However, this may be complicated by the possible presence of extragalactic 
magnetic fields, which would deflect the particles during their travel. Furthermore, since the GZK-energy losses are 
of stochastic nature, even a detailed knowledge of the extragalactic magnetic fields would not necessarily allow to 
follow a UHECR trajectory backwards to its source since the energy and therefor the Larmor radius of the particles 
have changed in an unknown way. Therefore it is not clear if charged particle astronomy with UHECRs is possible in 
principle or not. And even if possible, it remains unclear to which degree the angular resolution would be Umited by 
magnetic deflection. This topic will be discussed in Sect. 9.6. 



9.2. Severe Constraints on Scenarios producing more photons than hadrons 

The physics and astrophysics of UHECRs are also intimately Unked with the emerging field of neutrino astronomy 
(for reviews see Refs. [73]) as well as with the already established field of 7— ray astronomy (for reviews see, e.g., 
Ref. [74]). Indeed, all scenarios of UHECR origin, including the top-down models, are severely constrained by neutrino 
and 7— ray observations and limits. In turn, this Unkage has important consequences for theoretical predictions of fluxes 
of extragalactic neutrinos above about a TeV whose detection is a major goal of next-generation neutrino telescopes: If 
these neutrinos are produced as secondaries of protons accelerated in astrophysical sources and if these protons are not 
absorbed in the sources, but rather contribute to the UHECR flux observed, then the energy content in the neutrino flux 
can not be higher than the one in UHECRs, leading to the so called Waxman-Bahcall bound for transparent sources 
with soft acceleration spectra [20, 75]. This bound is shown in Fig. 2. If one of these assumptions does not apply, 
such as for acceleration sources with injection spectra harder than E^^ and/or opaque to nucleons, or in the top-down 
scenarios where X particle decays produce much fewer nucleons than y— rays and neutrinos, the Waxman-Bahcall 
bound does not apply, but the neutrino flux is still constrained by the observed diffuse 7— ray flux in the GeV range 
which is marked "EGRET" in Figs. 2 and 9. This bound whose implications will be discussed in the following section 
is marked "7— ray bound" in Figs. 2 and 3. 

Electromagnetic (EM) energy injected above the threshold for pair production on the CMB at ~ 10^^/ (1 eV at 
redshift z (to a lesser extent also on the infrared/optical background, with lower threshold) leads to an EM cascade, 
an interplay between pair production followed by inverse Compton scattering of the produced electrons. This cascade 
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continues until the photons fall below the pair production threshold at which point the universe becomes transparent for 
them. In todays universe this happens within just a few Mpc for injection up to the highest energies above lO^^eV. All 
EM energy injected above ~ 10^^ eV and at distances beyond a few Mpc today is therefore recycled to lower energies 
where it gives rise to a characteristic cascade spectrum c< E"^.! (jown to fractions of a GeV [76]. The universe thus 
acts as a calorimeter where the total EM energy injected above ^ lO'^/ (1 +z)eV is measured as a diffuse isotropic 
y— ray flux in the GeV regime. This diffuse flux is not very sensitive to the somewhat uncertain infrared/optical 
background [77]. Any observed diffuse y— ray background acts as an upper limit on the total EM injection. Since 
in any scenario involving pion production the EM energy fluence is comparable to the neutrino energy fluence, the 
constraint on EM energy injection also constrains allowed neutrino fluxes. 

This diffuse extragalactic GeV 7— ray background can be extracted from the total 7— ray flux measured by EGRET 
by subtracting the Galactic contribution. Since publication of the original EGRET limit in 1995 [78], models for 
this high latitude Galactic 7— ray foreground were improved significantly. This allowed the authors of Ref. [19] to 
reanalyze limits on the diffuse extragalactic background in the region 30 MeV-10 GeV and to lower it by a factor 
1.5-1.8 in the region around 1 GeV. There are even lower estimates of the extragalactic diffuse 7— ray flux [79]. In this 
article, however, we will use the more conservative limits from Ref. [19]. 

The energy in the extra-galactic 7— ray background estimated in Ref. [19] is slightly more than one hundred times 
the energy in UHECR above the GZK cutoff. The range of such trans-GZK cosmic rays is about ~ 30 Mpc, roughly 
one hundredth the Hubble radius, and only sources within that GZK range contribute to the trans-GZK cosmic rays. 
Therefore, any mechanism involving sources distributed roughly uniformly on scales of the GZK energy loss length 
~ 30Mpc and producing a comparable amount of energy in trans-GZK cosmic rays and photons above the pair 
production threshold can potentially explain this energy flux ratio. The details depend on the exact redshift dependence 
of source activity and other parameters and in general have to be verified by numerically solving the relevant transport 
equations, see, e.g., Ref. [6]. Such mechanisms include shock acceleration in powerful objects such as active galactic 
nuclei [80]. 

On the other hand, any mechanism producing considerably more energy in the EM channel above the pair production 
threshold than in trans-GZK cosmic rays tend to predict a ratio of the diffuse GeV 7— ray flux to the trans-GZK cosmic 
ray flux too high to explain both fluxes at the same time. As a consequence, if normalized at or below the observational 
GeV 7— ray background, such scenarios tend to explain at most a fraction of the observed trans-GZK cosmic ray 
flux. Such scenarios include particle physics mechanisms involving pion production by quark fragmentation, e.g. 
extra-galactic top-down mechanisms where UHECRs are produced by fragmenting quarks resulting from decay of 
superheavy relics [81]. Most of these quarks would fragment into pions rather than nucleons such that more 7— rays 
(and neutrinos) than cosmic rays are produced. Overproduction of GeV 7— rays can be avoided by assuming the sources 
in an extended Galactic halo with a high > 10^ overdensity compared to the average cosmological source density, 
which would also avoid the GZK cutoff [82]. These scenarios, however, start to be constrained by the anisotropy 
they predict because of the asymmetric position of the Sun in the Galactic halo for which there are no indications 
in present data [83]. Scenarios based on quark fragmentation also become problematic in view of a possible heavy 
nucleus component and of upper limits on the photon fraction of the UHECR flux [61]. 

As a specific example for scenarios involving quark fragmentation, we consider here the case of decaying Z- 
bosons. In this "Z-burst mechanism" Z-bosons are produced by UHE neutrinos interacting with the reUc neutrino 
background [84]. If the relic neutrinos have a mass my, Z-bosons can be resonantly produced by UHE neutrinos of 
energy Ey — M|/(2mv) — 4.2 x 10^' eV (eV/my). The required neutrino beams could be produced as secondaries of 
protons accelerated in high-redshift sources. The fluxes predicted in these scenarios have recently been discussed in 
detail, for example, in Refs. [85, 6]. In Fig. 9 we show an optimistic example taken from Ref. [6]. It is assumed that the 
relic neutrino background has no significant local overdensity. Furthermore, the sources are assumed to not emit any 
7— rays, otherwise the Z-burst model with acceleration sources over-produces the diffuse GeV 7— ray background [86]. 
We note that no known astrophysical accelerator exists that meets the requirements of the Z-burst model [86, 87]. 

However, a combination of new constraints discussed in the previous sections allows to rule out that the Z-burst 
mechanism explains a dominant fraction of the observed UHECR flux, even for pure neutrino emitting sources: 
As discussed in Sect. 8.4, a combination of cosmological data including the WMAP experiment limit the sum of 
the masses of active neutrinos to < leV [50]. In Sects. 8.1, 8.2 we have seen that solar and atmospheric neutrino 
oscillations indicate that individual neutrino masses are nearly degenerate on this scale [41], and thus the neutrino mass 
per flavor must satisfy niv < 0.33 eV. However, for such masses phase space constraints limit the possible over-density 
of neutrinos in our Local Group of galaxies to < 10 on a length scale of ~ 1 Mpc [88]. Since this is considerably smaller 
than the relevant UHECR loss lengths, neutrino clustering will not significantly reduce the necessary UHE neutrino 
flux compared to the case of no clustering. For the maximal possible value of the neutrino mass Wy — 0.33 eV, the 



Lectures on Astroparticle Physics 



February 2, 2008 



29 



CO 



E 
o 

> 

CM 

LU 

llT 



10^ 



10^ 



10 



EGRET 



10^ - 



-2 



1 Tv-r 



atm V 



1 \ \ \ \ \ 
GLUED 

~ORTE 



10= 




E[eV] 



FIGURE 9. From Ref. [6]. Flux predictions for a Z-burst model averaged over flavors and characterized by a neutrino injection 
flux per comoving volume oc E^^ up to 3 x 10^^ eV and for redshifts between and 3. The sources are assumed to be exclusive 
neutrino emitters. All neutrino masses were assumed equal with my = 0.33 eV and we again assumed maximal mixing between all 
flavors. The data and upper limits are as in Fig. 2. 



neutrino flux required for the Z-burst mechanism to explain the UHECR flux is only in marginal conflict with the 
FORTE upper limit [16], and factor 2 higher than the new GLUE limit [15], as shown in Fig. 9. For all other cases the 
conflict with both the GLUE and FORTE limits is considerably more severe. Also note that this argument does not 
depend on the shape of the low energy tail of the primary neutrino spectrum which could thus be even mono-energetic, 
as could occur in exclusive tree level decays of superheavy particles into neutrinos [89]. However, in addition this 
possibility has been ruled out by overproduction of GeV 7— rays due to loop effects in these particle decays [90]. 

The possibility that the observed UHECR flux is explained by the Z burst scenario involving normal astrophysical 
sources which produce both neutrinos and photons by pion production is already ruled out by the former EGRET limit: 
In this case the GeV 7— ray flux level would have roughly the height of the peak of the neutrino flux multiplied with 
the squared energy in Fig. 9, thus a factor ^ 100 higher than the EGRET level. 

Any further reduction in the estimated contribution of the true diffuse extra-galactic 7— ray background to the 
observed flux, therefore, leads to more severe constraints on the total EM injection. For example, future 7— ray 
detectors such as GLAST [91] will test whether the diffuse extragalactic GeV 7— ray background is truly diffuse 
or partly consists of discrete sources that could not be resolved by EGRET. Astrophysical discrete contributions such 
as from intergalactic shocks are in fact expected [92]. This could further improve the cascade limit to the point where 
even acceleration scenarios may become seriously constrained. 



9.3. New Primary Particles 

A possible way around the problem of missing counterparts within acceleration scenarios is to propose primary 
particles whose range is not limited by interactions with the CMB. Within the Standard Model the only candidate is 
the neutrino, whereas in extensions of the Standard Model one could think of new neutrals such as axions or stable 
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supersymmetric elementary particles. Such options are mostly ruled out by the tension between the necessity of a 
small EM coupling to avoid the GZK cutoff and a large hadronic coupling to ensure normal air showers [93]. Also 
suggested have been new neutral hadronic bound states of light gluinos with quarks and gluons, so-called R-hadrons 
that are heavier than nucleons, and therefore have a higher GZK threshold [94], as can be seen from Eq. (78). Since 
this too seems to be disfavored by accelerator constraints [95] we will here focus on neutrinos. 

In both the neutrino and new neutral stable particle scenario the particle propagating over extragalactic distances 
would have to be produced as a secondary in interactions of a primary proton that is accelerated in a powerful active 
galactic nucleus which can, in contrast to the case of EAS induced by nucleons, nuclei, or 7— rays, be located at 
high redshift. Consequently, these scenarios predict a correlation between primary arrival directions and high redshift 
sources. In fact, possible evidence for a correlation of UHECR arrival directions with compact radio quasars and BL- 
Lac objects, some of them possibly too far away to be consistent with the GZK effect, was recently reported [96]. 
The main challenge in these correlation studies is the choice of physically meaningful source selection criteria and the 
avoidance of a posteriori statistical effects. However, a moderate increase in the observed number of events will most 
likely confirm or rule out the correlation hypothesis. Note, however, that these scenarios require the primary proton to 
be accelerated up to at least 10^^ eV, demanding a very powerful astrophysical accelerator. 



9.4. New Neutrino Interactions 

Neutrino primaries have the advantage of being well established particles. However, within the Standard Model 
their interaction cross section with nucleons shown in Fig. 1 falls short by about five orders of magnitude to produce 
air showers starting high in the atmosphere as observed. Electroweak instantons could change this but this possibility 
is speculative [97]. The neutrino-nucleon cross section, cTva?, however, can be enhanced by new physics beyond the 
electroweak scale in the CM frame, or above about a PeV in the nucleon rest frame. Note that the CM energy reached 
by an UHECR nucleon of energy E interacting with an atmospheric nucleon at rest is ^/s ~ 0.4(£'/10^'^eV)PeV. 
Neutrino induced air showers may therefore rather directly probe new physics beyond the electroweak scale. 

One possibility consists of a large increase in the number of degrees of freedom above the electroweak scale [98]. 
A specific instance of this idea appears in theories with n additional large compact dimensions and a quantum gravity 
scale M4_|_„ ~ TeV that has recently received much attention in the literature [99] because it provides an alternative 
solution to the hierarchy problem in grand unifications of gauge interactions without a need of supersymmetry. The 
idea is to dimensionally reduce the n + 4 dimensional gravitational action 

(123) 

with g the determinant of the metric and R the Ricci scalar, to four dimensions by integrating out the n compact 
dimensions. This yields the relation 

Ml, = VnMlXl, (124) 

where Mpi and Vn are the four-dimensional Planck mass and the volume of the n extra dimensions, respectively. The 
weakness of gravity can now be understood as a consequence of the fact that it is the only force that propagates into the 
extra dimensions: Their large volume dilutes gravitational interactions between the Standard Model particles which 
are confined to a 3-brane representing our world. For compact extra dimensions, gravity would only be modified at 
scales below 

/ Mpi \ 2/" r TeV \ / Mp, \ 




r„ ~ M7, „ — 10"" — — cm, (125) 

which is r„ < 1mm for n > 2, M4+„ >TeV, and thus consistent with gravity tests at small distances. In contrast, 
non-gravitational interactions are confined to the 3-brane and thus the Standard Model is not modified. 

The neutrino-nucleon cross section in these frameworks is obtained by substituting the new fundamental cross 
section C7,(xi, Q) for the electroweak cross section at the parton level in Eq. (75), 

= '£fi{x,Q)Oi{xs,Q) . (126) 

One of the largest contributions to the neutrino-nucleon cross section turns out to be the production on our 3-brane of 
microscopic black holes which are solutions of 4 -|- n-dimensional gravity described by Eq. (123). These cross sections 
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scale as {s/ Mj^^ for s > Mj^^ . Their UV-divergence is due to the non-renormalizable, classical character of 
the gravitational interaction Eq. (123) in the sense of Sect. 3. The production of compact branes, completely wrapped 
around the extra dimensions, may provide even larger contributions [100]. The resulting total cross sections can be 
larger than in the Standard Model by up to a factor ~ 100 if M4+„ ~TeV [101]. However, extra dimensions with a 
flat geometry are severely constrained by astrophysics: Core collapse of massive stars would lead to production of 
gravitational excitations in the large compact extra dimensions, mostly by nucleon-nucleon bremsstrahlung. These so 
called Kaluza- Klein gravitons of mass nig are then gravitationally trapped around the newly bom neutron star during 
their livetime T ^Mpj/m^, ^ 10'^(10MeV/;7i^,)^ yr Their subsequent decay into two y— rays would make neutron stars 
shine in 7— rays. The non-observation of such emission leads to lower bounds on M4+„ which decrease with increasing 
n, starting with M5 > lO^TeV and going down to Mg > 1 TeV and still lower values for larger n [102, 103], for flat 
compact extra dimensions. This implies that significant contributions to the neutrino-nucleon cross section in these 
extra dimension scenarios require either n>5 extra dimensions or a warped geometry. 

Whereas the sub-hadronic scale cross sections obtained in some extra dimension scenarios are still too small to 
be consistent with observed air showers and thus to explain the observed UHECR events [104], they can stiU have 
important phenomenological consequences. This is because UHECR data can be used to put constraints on cross 
sections satisfying Ovt^{E > lO'^eV) < 10^^^ cm^. Particles with such cross sections would give rise to horizontal air 
showers which have not yet been observed. Resulting upper limits on their fluxes assuming the Standard Model cross 
section Eq. (7) are shown in Fig. 2. Comparison with the "cosmogenic" neutrino flux produced by UHECRs interacting 
with the CMB then results in upper limits on the cross section which are about a factor 1000 larger than Eq. (7) in 
the energy range between ~ lO'^eV and ~ lO^'^eV [105, 106, 107]. The projected sensitivity of future experiments 
shown in Fig. 3 indicate that these limits could be lowered down to the Standard Model cross section [107]. In case of 
a detection of penetrating events the degeneracy of the cross section with the unknown neutrino flux could be broken 
by comparing the rates of horizontal air showers with the ones of Earth skimming events [108]. This would allow to 
"measure" the neutrino-nucleon cross section at energies unreachable by any forseeable terrestrial accelerator ! 



9.5. Violation of Lorentz Invariance 

The most elegant solution to the problem of apparently missing nearby sources of UHECRs and for their putative 
correlation with high redshift sources would be to speculate that the GZK effect does not exist theoretically. A number 
of authors pointed out [109, 110] that this may be possible by allowing violation of Lorentz invariance (VLl) by a tiny 
amount that is consistent with all current experiments. At a purely theoretical level, several quantum gravity models 
including some based on string theories do in fact predict non-trivial modifications of space-time symmetries that also 
imply VLI at extremely short distances (or equivalently at extremely high energies); see e.g., Ref. [Ill] and references 
therein. These theories are, however, not yet in forms definite enough to allow precise quantitative predictions of 
the exact form of the possible VLl. Current formulations of the effects of a possible VLl on high energy particle 
interactions relevant in the context of UHECR, therefore, adopt a phenomenological approach in which the form of 
the possible VLI is parametrized in various ways. VLI generally implies the existence of a universal preferred frame 
which is usually identified with the frame that is comoving with the expansion of the Universe, in which the CMB is 
isotropic. 

A direct way of introducing VLI is through a modification of the standard dispersion relation, — = m^, 
between energy E and momentum p = \p\ of particles, m being the invariant mass of the particle. Currently there is no 
unique way of parameterizing the possible modification of this relation in a Lorentz non-invariant theory. We discuss 
here a parameterization of the modified dispersion relation which covers most of the qualitative cases discussed in the 
Uterature and, for certain parameter values, aUows to completely evade the GZK Umit, 

E^-p2-m^^-2dE^-^—-i;K- (127) 

Here, the Planck mass Mpi characterizes non-renormalizable effects with dimensionless coefficients ^ and ^, and the 
dimensionless constant d exemplifies VLI effects due to renormalizable terms in the Lagrangian, see the discussion in 
Sect. 3. The standard Lorentz invariant dispersion relation is recovered in the limit £,,^,d ^ 0. 

The constants dj^O can break Lorentz invariance spontaneously when certain Lorentz tensors c^v have coupUngs 
to fermions of the form d^vW7'^^^¥^ and acquire vacuum expectation values of the form {d^v) = d5^5y. If rota- 
tional invariance and gauge symmetry are preserved, such renormalizable Lorentz invariance breaking terms in the 
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Lagrangian, whose Lorentz invariant part is given by Eq. (58), are characterized by a single time-like vector u^, with 
u^Ufi = — 1 , which defines a preferred reference frame [112]. The dimensionless terms can be interpreted as a change 
of the maximal particle velocity [110] Vmax = dE/dp\E,p:»m — 1—d. At a fixed energy £ one has the correspondence 
d^{^ /2) {E /Mpi) + (C /2) {E /Mpi)2, as can be seen from Eq. (127). 

Within effective field theory, effects of first order in ^ 7^ 0, arise from the most general terms of the form 

2^«%v(« • d)u^F^'' + ^u^xjfY^(?,,+X275){u ■ 5) V, (128) 

where F denotes the dual of the field strength F. For photons and electrons this leads to E, = ±K and £, = Xi± 
hi, respectively, in Eq. (127), where ± refers to helicity which remains conserved in the presence of the terms 
Eq. (128) [113]. These terms also violate CPT . Effects of first order in Mpj^ 7^ 0, are possible, for example, in 
non-critical Liouville string theory due to recoiling D-branes [114]. Finally, in critical string theory, effects second 
order in M^^ , C 7^ 0, can be induced due to quantum gravity effects. 

Interestingly, it has been pointed out recently that in the supersymmetric Standard Model VLI terms must be non- 
renormalizable and do not lead to modifications of any dispersion relations [115]. 

Now, consider the GZK photo-pion production process in which a nucleon of energy E, momentum p and mass 
niN collides head-on with a CMB photon of energy e producing a pion and a recoiling nucleon. The threshold initial 
momentum of the nucleon for this process according to standard Lorentz invariant kinematics is 

Pth,o = (w^ -|-2m;tmAf)/4e , (129) 

where and m^v are the pion and nucleon masses, respectively. Assuming exact energy-momentum conservation but 
using the modified dispersion relation given above, in the ultra-relativistic regime m <C f <C M, and neglecting sub- 
leading terms, the new nucleon threshold momentum />th under the modified dispersion relation Eq. (127) for J = 
satisfies [116] 

-/3/-ax'-|-x-l=0, (130) 

where x = pth/ Pth.o. and 

^ ^ ^^^fe,0 mj^mN_^ ^^^^^ 



nijiinN 



2{ml + 2mnmN)M^ {m^+rriNY ' 

One can show that the same modified dispersion relation Eq. (127) leads to the same condition Eq. (130) for absorption 
of high energy gamma rays through e+e~ pair production on the infrared, microwave or radio backgrounds, if one 
substitutes pth,o = m^/e, a = ^ Ptho/ (S'WeAfpi), P = ^Cfth o/(16»J^Mp[), where ntg is the electron mass. 

If 1^ , ~ 1, there is no real positive solution of Eq. (130), implying that the GZK process does not take place and 
consequently the GZK cutoff effect disappears completely. Thus UHE nucleons and/or photons will be able to reach 
Earth from any distance. On the other hand, if future UHECR data confirm the presence of a GZK cutoff at some 
energy then that would imply upper limits on the couplings ^ and ^, thus probing specific Lorentz non-invariant 
theories. If pth — Pth.o, one could conclude from Eq. (130) that a,j3 < 1, which translates into \^ \ < 10"'^ for the first 
order effects, and | ^ | < 10~^ for the second order effects, ^ = [1 16]. These values correspond to values \d\< 10~^^ 
for the paremeters of renormalizable VLI. Confirmation of a cut-off for TeV photons with next-generation 7— ray 
observatories would lead to somewhat weaker constraints [117]. 

More generally, modification of reaction kinematics or new reaction channels are expected whenever the terms on 
the right hand side of Eq. (127) become comparable to m^, in rough numbers. 
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Note that by far the smallest parameter values would be probed by particles with the smallest mass, specifically the 
neutrino, m <eV at the highest energies. This makes the prospects of future detections of cosmogenic neutrinos, see, 
e.g.. Fig. 2, very exciting also for VLI constraints. 

In addition, the non-renormaUzable terms in the dispersion relation Eq. (127) imply a change in the group velocity 
which for the first-order term leads to time delays over distances r given by 

A;~fr— ) ( — I sec. (133) 

For 1^ I ~ 1 such time delays could be measurable, for example, by fitting the arrival times of 7— rays arriving from 
7— ray bursts to the predicted energy dependence. 

We mention that if VLI is due to modification of the space-time structure expected in some theories of quantum 
gravity, for example, then the strict energy-momentum conservation assumed in the above discussion, which requires 
space-time translation invariance, is not guaranteed in general, and then the calculation of the modified particle 
interaction thresholds becomes highly non-trivial and non-obvious. Also, it is possible that a Lorentz non-invariant 
theory while giving a modified dispersion relation also imposes additional kinematic structures such as a modified 
law of addition of momenta. Indeed, Ref. [Ill] gives an example of a so-called K"-Minkowski non-commutative 
space- time in which the modified dispersion relation has the same form as in Eq. (127) but there is also a modified 
momentum addition rule which compensates for the effect of the modified dispersion relation on the particle interaction 
thresholds discussed above leaving the threshold momentum unaffected and consequently the GZK problem unsolved. 
In scenarios where the relativity of inertial frames is preserved by a non-linear representation of the Poincare group, 
thresholds are in general significantly modified only if the effective mass scale Mpi/ ^ is of the order of the unmodified 
threshold energy in the laboratory frame [118]. 

There are several other fascinating effects of allowing a small VLI, some of which are relevant for the question 
of origin and propagation of UHECR. For example, any movement relative to the preferred frame defined by 
in Eq. (128) gives rise to spatial anisotropy. Clock comparison and spin precession experiments then lead to limits 
on the dimensionless parameters in Eq. (128) between and i^(10^^), depending on the particle [113]. Similar 
Umits result from astrophysical arguments: The observation of polarized MeV synchrotron radiation from electrons 
in the Crab nebula implies the absence of vacuum Cerenkov radiation e — > ey for electrons up to energies E ~ 
l.SPeV [119, 120]. This process can become possible if the electron speed becomes larger than the speed of light at 
high energies and leads to limits on VLI of size comparable to the before mentioned laboratory constraints [119, 120]. 
These constraints basically rule out effects of order £/Mpi which might be a challenge for certain quantum gravity 
scenarios [121]. Note that these current constraints on VLI parameters still allow strong modification of GZK 
kinematics by VLI parameters of the order given in Eq. (132). 



9.6. Cosmic IMagnetic Fields and Their Influence on Ultra-High Energy Cosmic Ray 

Propagation 

Cosmic magnetic fields are inextricably linked with cosmic rays in several respects. First, they play a central role in 
Fermi shock acceleration. Second, large scale extra-galactic magnetic fields (EGMF) can cause significant deflection 
of charged cosmic rays during propagation and thus obviously compUcate the relation between observed UHECR 
distributions and their sources. 

Magnetic fields are omnipresent in the Universe, but their true origin is still unclear [122]. Magnetic fields in galaxies 
are observed with typical strengths of a few micro Gauss, but there are also some indications for fields correlated 
with larger structures such as galaxy clusters [123]. Magnetic fields as strong as ~ 1/iG in sheets and filaments of 
the large scale galaxy distribution, such as in our Local Supercluster, are compatible with existing upper limits on 
Faraday rotation [123, 124, 125]. It is also possible that fossil cocoons of former radio galaxies, so called radio ghosts, 
contribute significantly to the isotropization of UHECR arrival directions [126]. 

To get an impression of typical deflection angles one can characterize the EGMF by its r.m.s. strength B and a 
coherence length Ic. If we neglect energy loss processes for the moment, then the r.m.s. deflection angle over a distance 
r > /c in such a field is (£,?•) ~ (2r/(;/9)'/^/rL [127], where the Larmor radius of a particle of charge Ze and energy 
£■ is ~ E/{ZeB). In numbers this reads 

^(^.^)-0.8°z(^) '(t^)"'(t^)"'(t^), (134) 
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for r > Ic. This expression makes it immediately obvious that fields of fractions of micro Gauss lead to 
strong deflection even at the highest energies. This goes along with a time delay x{E,r) ~ rQ{E,d)'^ ~ 
1.5 X 102z2(£:/1020eV)-2(r/10Mpc)2(;c/Mpc)(fi/10-''G)2yr which can be milHons of years. A source visi- 
ble in UHECRs today could therefore be optically invisible since many models involving, for example, active galaxies 
as UHECR accelerators, predict variability on shorter time scales. 

Quite a few simulations of the effect of extragalactic magnetic fields (EGMF) on UHECRs exist in the hterature, 
but usually idealizing assumptions concerning properties and distributions of sources or EGMF or both are made: 
In Refs. [128, 129, 130, 131, 132] sources and EGMF follow a pancake profile mimicking the local supergalactic 
plane. In other studies EGMF have been approximated in a number of fashions: as neghgible [133, 134], as stochastic 
with uniform statistical properties [135, 136, 137], or as organized in spatial cells with a given coherence length 
and a strength depending as a power law on the local density [138]. Only recently attempts have been made to 
simulate UHECR propagation in a realistically structured universe [139, 140]. For now, these simulations are limited 
to nucleons. 

In Ref. [139] the magnetized extragalactic environment used for UHECR propagation is produced by a simulation 
of the large scale structure of the Universe. The simulation was carried out within a computational box of 50/z^' Mpc 
length on a side, with normalized Hubble constant /z = //o/(100 km s^' Mpc^' ) = 0.67, and using a comoving grid 
of 512^ zones and 256^ dark matter particles. The EGMF was initiahzed to zero at simulation start and subsequently 
its seeds were generated at cosmic shocks through the Biermann battery mechanism [141]. Since cosmic shocks form 
primarily around collapsing structures including filaments, the above approach avoids generating EGMF in cosmic 
voids. 

In Ref. [140] constrained simulations of the local large scale structure were performed and the magnetic smoothed 
particle hydrodynamics technique was used to follow EGMF evolution. The EGMF was seeded by a uniform seed 
field of maximal strength compatible with observed rotation measures in galaxy clusters. 

The questions considered in these two works were somewhat different, however. In Ref. [140] deflections of 
UHECR above 4 x 10^'^eV were computed as a function of the direction to their source which were assumed to 
be at cosmological distances. This made sense, because (i) the constrained simulations gives a viable model of our 
local cosmic neighborhood within about 100 Mpc, at least on scales beyond a few Mpc and (ii) the deflections typically 
were found to be smaller than a few degrees. Concrete source distributions were not considered. 

In contrast, Ref. [139] was not concerned with concrete sky distributions or deflection maps because the simulation 
was unconstrained and thus only gave a typical large scale structure model and not our concrete local neighborhood. 
Instead, the question was asked which observer positions and source distributions and characteristics lead to UHECR 
distributions whose spherical multi-poles for / < 10 and auto-correlation at angles < 20° are consistent with 
observations. As a result it was found that (i) the observed large scale UHECR isotropy requires the neighborhood 
within a few Mpc of the observer is characterized by weak magnetic fields below 0.1 /iG, and (ii) once that choice 
is made, current data do not strongly discriminate between uniform and structured source distributions and between 
negligible and considerable deflection. Nevertheless, current data moderately favor a scenario in which (iii) UHECR 
sources have a density Us ~ lO^^Mpc"^ and follow the matter distribution and (iv) magnetic fields are relatively 
pervasive within the large scale structure, including filaments, and with a strength of order of a \iQ in galaxy clusters. 
A two-dimensional cut through the EGMF environment of the observer in a typical such scenario is shown in Fig. 10. 

It was also studied in Ref. [139] how future data of considerably increased statistics can be used to learn more 
about EGMF and source characteristics. In particular, low auto-correlations at degree scales imply magnetized sources 
quite independent of other source characteristics such as their density. The latter can only be estimated from the auto- 
correlations halfway reliably if magnetic fields have negligible impact on propagation. This is because if sources are 
immersed in considerable magnetic fields, their images are smeared out, which also smears out the auto-correlation 
function over several degrees. For a sufficiently high source density, individual images can thus overlap and sensitivity 
to source density is consequently lost. The statistics expected from next generation experiments such as Pierre 
Auger [72] and EUSO 1 142| should be sufficient to test source magnetization by the auto-correlation function [139]. 

Interestingly, however, there is a considerable quantifiable difference in the typical deflection angles predicted by 
the two EGMF scenarios in Refs. [139, 140] that can not be compensated by specific source distributions: Even for 
homogeneous source distributions, the average deflection angle for UHECRs above 4 x lO'^eV obtained in Ref. [139] 
is much larger than in Ref. [140], as can be seen in Fig. 11. In fact, even if the magnetic field strength is reduced 
by a factor 10 in the simulations of Ref. [139], the average deflection angle above 4 x 10^^ eV is stiU ~ 30°, only 
a factor 2.2 smaller. This non-linear behavior of deflection with field normaUzation is mostly due to the strongly 
non-homogeneous character of the EGMF. 
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FIGURE 10. Log-scale two-dimensional cut through magnetic field total strength in Gauss (color scale in Gauss) for a scenario 
in good agreement with UHECR data studied in Ref. [139]. The observer is in the center of the figures and is marked by a star. The 
EGMF strength at the observer is ~ 10^" G. 



Most of these differences are probably due to the different numerical models for the magnetic fields. Although 
Ref. [140] start with uniform seed fields, whereas in Ref. [139] seed fields are injected at shocks, by itself, this 
difference should not influence the resulting EGMF very much at late times, at least inside galaxy clusters [124]. It 
should be noted, however, that in the filaments, where the gas motions are more uniform, the simulated magnetic 
fields may depend to a certain extent on the initial seed fields although that is not trivial to quantify in general terms. 
In addition, numerical resolution may play an important role because it affects the amplification and the topological 
structure of the magnetic fields, both of which are important for the normalization procedure, see below. The resolution 
in Ref. [139] is constant and much better in filaments and voids but worse in the core of galaxy clusters than the 
(variable) resolution in Ref. [140]. If in both simulations the magnetic fields are normalized to (or reproduce) the same 
"observed" values in the core of rich clusters then obviously their values in the filaments will be very different for 
the reasons outlined above. This may partly explain why the contribution of filaments to UHECR deflection is more 
important in Ref. [139], although a more detailed analysis and comparison are required to settle the issue. In any case, 
the magnetic fields obtained in Ref. [139] seem to be quite extended, as can be seen in Fig. 12: About 10% of the 
volume is filled with fields stronger than 10 nano Gauss, and a fraction of 10^^ is filled by fields above a micro Gauss. 
Furthermore, typical deflection angles change at most by a factor of 2 if magnetic field normalization is decreased by 
a factor 10 or seed fields are chosen as uniform in these simulations. The different amounts of deflection obtained in 
the simulations of Refs. [139, 140] show that the distribution of EGMF and their effects on UHECR propagation are 
currently rather uncertain. 

Finally we note that these studies should be extended to include heavy nuclei [143] since there are indications 
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FIGURE 11. The cumulative distribution of UHECR deflection angles a with respect to the line of sight to the sources. This 
is for a scenario from Ref. [139] in good agreement with UHECR data, where the sources follow the baryon density and have 
average density tis — 2.4 x lO^^Mpc^-', and the EGMF included in the large scale structure simulation reaches several micro 
Gauss in the most prominent galaxy cluster. Shown are the average (middle, histogram) and 1-ct variations (upper and lower 
curves) above 4 x lO'** eV, over 24 realizations varying in the positions and luminosities Qi of individual sources, the latter assumed 
to be distributed as drisjclQi = QJ^'^ with 1 < 2, < 100 in arbitrary units. Also given on top of the figure are average and variances 
of the distributions. 



that a fraction as large as 80% of iron nuclei may exist above 10^ eV [61]. As a consequence, even in the EGMF 
scenario of Ref. [140] deflections could be considerable and may not allow particle astronomy along many lines of 
sight: The distribution of deflection angles in Ref. [140] shows that deflections of protons above 4 x lO'^eV of > 1° 
cover a considerable fraction of the sky. Suppression of deflection along typical Unes of sight by small filling factors 
of deflectors is thus unimportant in this case. The deflection angle of any nucleus at a given energy passing through 
such areas will therefore be roughly proportional to its charge as long as energy loss lengths are larger than a few 
tens of IVIpc [144]. Deflection angles of ~ 20° at ~ 4 x 10^^ eV should thus be the rule for iron nuclei. In contrast 
to the contribution of our Galaxy to deflection which can be of comparable size but may be corrected for within 
sufficiently detailed models of the galactic field, the extra-galactic contribution would be stochastic. Statistical methods 
are therefore hkely to be necessary to learn about UHECR source distributions and characteristics. In addition, should a 
substantial heavy composition be experimentally confirmed up to the highest energies, some sources would have to be 
surprisingly nearby, within a few Mpc, otherwise only low mass spallation products would survive propagation [145]. 

The putative clustered component of the UHECR flux whose fraction of the total flux seems to increase with 
energy [65] may play a key role in this context. It could be caused by discrete sources in directions with small 
deflection. Spectrum and composition of the flux from such sources could still by modified considerably by magnetic 
fields concentrated around the source [131, 146]. For example, since, apart from energy losses, cosmic rays of same 
rigidity Z/A are deflected similarly by cosmic magnetic fields, one may expect that the composition of the clustered 
component may become heavier with increasing energy. Indeed, in Ref. [147] it was speculated that the AGASA 
clusters may be consistent with consecutive He, Be-Mg, and Fe bumps. 

It thus seems evident that the influence of large scale cosmic magnetic fields on ultra-high energy cosmic ray 
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FIGURE 12. The cumulative filling factors for EGMF strength in the simulations used in Ref. [139] above (decreasing curve) 
and below (increasing curve) a given threshold, as a function of that threshold. 

propagation is currently hard to quantify and may not allow to do "particle astronomy" along most lines of sight, 
especially if a significant heavy nucleus component is present above lO'^ eV. 

ACKNOWLEDGMENTS 

The material presented here is based on a university course on neutrino physics taught by the author and on research 
work with various collaborators of whom I would especially hke to thank Torsten Ensshn, Francesco Miniati and 
Dmitry Semikoz. Finally, I would hke to thank the organizers of the Xlth Brazilian School of Cosmology and 
Gravitation for a terrific school. 



Lectures on Astroparticle Physics 



February 2, 2008 



38 



REFERENCES 



1. D. H. Perkins (1987). Introduction to High Energy Physics, Addison- Wesley. 

2. S. Weinberg (1972). Gravitation and Cosmology, New York: John Wiley & Sons. [A detailed exposition of general relativity 
and its applications]. 

3. S. Weinberg (1995). The Quantum Theory of Fields, volume 1: Foundations, Cambridge: Cambridge University Press. [Basic 
concepts of relativistic quantum field theory with special emphasize on how its general form is determined by symmetry 
principles: scattering theory, Feynman rules, quantum electrocl3Tiamics, path integrals, renormalization]. 

4. S. Weinberg (1996). The Quantum Theory of Fields, volume 2: Modern Applications, Cambridge: Cambridge University 
Press. [Introduction to non-abelian gauge theories: renormalization group, spontaneous symmetry breaking, anomalies, 
non-perturbative field configurations]. 

5. R. Gandhi, C. Quigg, M. H. Reno and I. Sarcevic, Phys. Rev. D 58, 093009 (1998) [arXiv:hep-ph/9807264]. 

6. D. V. Semikoz and G. Sigl, JCAP 0404, 003 (2004) [arXiv:hep-ph/0309328]. 

7. I acknowledge Dmitry Semikoz for updating this figure. 

8. see, e.g., http://ast.leeds.ac.uk/ workshop/halzen_leeds.ppt. 

9. For general information see http://wsgs02.lngs.infn. it:8000/macro/; see also M. Ambrosio et al. [MACRO Collaboration], 
Astropart. Phys. 19, 1 (2003) [arXiv:astro-ph/0203181]. 

10. P Niessen [the AMANDA Collaboration], arXiv:astro-ph/0306209; For the energy range 2.5 x lO'^ cV < £ < 5.6 x lO'^ eV 
we reseated the AMANDA-BIO limit from http://www-rccn.icrr.u-tokyo.ac.jp/icrc2003/PROCEEDINGS/PDF/326.pdf 
to AMANDA-II exposure. 

11. V. Balkanov et al. [BAIKAL Collaboration], Nucl. Phys. Proc. Suppl. 110, 504 (2002) [arXiv:astro-ph/0112446]; for general 
information see http://www-zeuthen.desy.de/baikal/baikalhome.html 

12. S. Yoshida for the AGASA Collaboration, Proc. of 27th ICRC (Hamburg) 3, 1 142 (2001). 

13. R. M. Baltrusaitis et al, Astrophys. J. 281, L9 (1984) ; Phys. Rev. D 31, 2192 (1985). 

14. I. Kravchenko et al, Astropart. Phys. 20, 195 (2003) [arXiv:astro-ph/0206371]; I. Kravchenko, arXiv:astro-ph/0306408. 

15. P. W. Gorham, K. M. Liewer and C. J. Naudet, astro-ph/9906504; P. W. Gorham, K. M. Liewer, C. J. Naudet, D. P. Saltzberg 
and D. R. Williams, arXiv:astro-ph/0102435; P W. Gorham, C. L. Hebert, K. M. Liewer, C. J. Naudet, D. Saltzberg and 
D. Williams, Phys. Rev. Lett. 93, 041101 (2004) [arXiv:astro-ph/0310232]. 

16. N. G. Lehtinen, R W. Gorham, A. R. Jacobson and R. A. Roussel-Dupre, Phys. Rev. D 69, 013008 (2004) [arXiv:astro- 
ph/0309656]. 

17. M. Takeda et al., Phys. Rev. Lett. 81, 1163 (1998) [arXiv:astro-ph/9807193]; M. Takeda et al, Astrophys. J. 522, 
225 (1999) [arXiv:astro-ph/9902239]; N. Hayashida et al, arXiv:astro-ph/0008102; see also http ://www-akeno.lcrr.u- 
tokyo.ac.jp/AGASA/. 

18. R. U. Abbasi et al. [High Resolution Fly's Eye Collaboration], Phys. Rev. Lett. 92, 151101 (2004) [arXiv:astro-ph/0208243]; 
T. Abu-Zayyad et al. [High Resolution Fly's Eye Collaboration], arXiv:astro-ph/0208301; D. R. Bergman [The HiRes 
Collaboration], arXiv:astro-ph/0407244. 

19. A. W. Strong, 1. V. Moskalenko and O. Reimer, arXiv:astro-ph/0306345. 

20. E. Waxman and J. N. Bahcall, Phys. Rev D 59, 023002 (1999) [arXiv:hep-ph/9807282]; J. N. Bahcall and E. Waxman, Phys. 
Rev. D 64, 023002 (2001) [arXiv:hep-ph/9902383]. 

21. J. J. Blanco-Pillado, R. A. Vazquez and E. Zas, Phys. Rev. Lett. 78, 3614 (1997) [astro-ph/9612010]; K. S. Capelle, 

J. W. Cronin, G. Parente and E. Zas, Astropart. Phys. 8, 321 (1998) [astro-ph/9801313]; A. Letessier-Selvon, Nucl. Phys. Proc. 
Suppl. 91, 473 (2000) [astro-ph/0009416]; X. Bertou, R Billoir, O. Deligny, C. Lachaud and A. Letessier-Selvon, Astropart. 
Phys. 17, 183 (2002) [astro-ph/0104452]. 

22. M. Sasaki and M. Jobashi, Astropart. Phys. 19, 37 (2003) [arXiv:astro-ph/0204167]. 

23. G. W. S. Hou and M. A. Huang, Nucl. Phys. Proc. Suppl. 118, 425 (2003) [arXiv:astro-ph/0204145]. 

24. see S. Bottai and S. Giurgola, http://www-rccn.icrr.u-tokyo.ac.jp/icrc2003/PROCEEDINGS/PDF/279.pdf. 

25. For general information see http://antares.in2p3.fr; see also ANTARES Collaboration, e-print astro-ph/9907432; T. Montaruli 
etal. [ANTARES Collaboration], arXiv:physics/0306057. 

26. For general information see http://www.nestor.org.gr . See also L. Resvanis, Proc. Int. Workshop on Neutrino Telescopes, 
Venice 1999, vol. II, 93. 

27. For general information see http://icecube.wisc.edu/; see also F. Halzen: Am. Astron. Soc. Meeting 192, # 62 28 (1998); 
AMANDA collaboration: astro-ph/9906205, Proc. 8''' International Workshop on Neutrino Telescopes, Venice, Feb. 1999. 

28. see, e.g., talk by U. Katz on HENA 2003 workshop, http://antares.in2p3.fr/antares/stolar/km3/program.htm. 

29. I. Kravchenko et al, Astropart. Phys. 20, 195 (2003) [arXiv:astro-ph/0206371]; for general information on RICE see 
http://kuhep4.phsx.ukans.edu/iceman/index.html. 

30. P. Gorham et al. (ANITA collaboration), http://www.ps.uci.edu/ barwick/an itaprop.pdf. 

31. R. N. Mohapatra, P. B. Pal (2004). Massive Neutrinos in Physics and Astrophysics, Imperial College Press. 

32. H. V. Klapdor-Kleingrothaus, A. Dietz, I. V. Krivosheina, C. Dorr and C. Tomei [Heidelberg-Moscow-Genius collaborations], 
Phys. Lett. B 578, 54 (2004) [arXiv:hep-ph/0312171]. 

33. see, e.g., C. Amaboldi et al. [CUORE Collaboration], Astropart. Phys. 20, 91 (2003) [arXiv:hep-ex/0302021]; C. Amaboldi 
et al. [CUORE Collaboration], Nucl. Instrum. Meth. A 518, 775 (2004) [arXiv:hep-ex/0212053]. 

34. C. Weinheimer, Nucl. Phys. Proc. Suppl. 118 (2003) 279. 

35. see http://www-ikl.fzk.de/tntiuw/. 



Lectures on Astroparticle Physics 



February 2, 2008 



39 



36. G. G. Raifelt (1996). Stars as Laboratories for Fundamental Physics, University of Chicago Press. 

37. J. N. Bahcall and M. H. Pinsonneault, Phys. Rev. Lett. 92, 121301 (2004) [arXiv:astro-ph/0402114]. 

38. see http://www.sno.phy.qucensu.ca/. 

39. see http://www.awa.tohoku.ac.jp/KamLAND/. 

40. J. N. Bahcall, M. C. Gonzalez-Garcia and C. Pena-Garay, arXiv:hep-ph/0406294. 

41. M. Maltoni, T. Schwetz, M. A. Tortola and J. W. F. Valle, arXiv:hep-ph/0405172. 

42. see http://neutrino.kek.jp/. 

43. see http://www-sk.icrr.u-tokyo.ac.jp/doc/sk/. 

44. Y. Ashie et al. [Super-Kamiokande Collaboration], arXiv:hep-ex/0404034. 

45. see, e.g., G. Altarelli and F. Feruglio, arXiv:hep-ph/0405048. 

46. F R. Klinkhamer, arXiv:hep-ph/0407200. 

47. E. W. Kolb and M. S. Turner (1990). The Early Universe, Addison- Wesley 

48. D. N. Spergel et al, Astrophys. J. Suppl. 148, 175 (2003) [arXiv:astro-ph/0302209]. 

49. G. Steigman, arXiv:hep-ph/0309347. 

50. S. Hannestad, arXiv:hep-ph/0404239; S. Hannestad, JCAP 0305, 004 (2003) [arXiv:astro-ph/0303076]. 

51. S. W. Allen, R. W. Schmidt and S. L. Bridle, Mon. Not. Roy Astron. Soc. 346, 593 (2003) [arXiv:astro-ph/0306386]. 

52. U. Seljak et al, "Cosmological parameter analysis including SDSS Ly-alpha forest and galaxy arXiv:astro-ph/0407372. 

53. A. D. Dolgov, S. H. Hansen, S. Pastor, S. T. Petcov, G. G. Raffelt and D. V. Semikoz, Nucl. Phys. B 632, 363 (2002) 
[arXiv:hep-ph/0201287]. 

54. W. Buchmuller, P. Di Bari and M. Plumacher, arXiv:hep-ph/0406014. 

55. A. D. Sakharov, Pisma Zh. Eksp. Teor. Fiz. 5, 32 (1967) [JETP Lett. 5, 24 (1967 SOPUA,34,392-393.1991 UFNAA,161,61- 
64.1991)]. 

56. V. F Hess, Phys. Z. 13 (1912) 1084. 

57. P Auger, R. Maze, T. Grivet-Meyer, Academie des Sciences 206 (1938) 1721; P Auger, R. Maze, ibid. 207 (1938) 228. 

58. for a general introduction on cosmic rays see, e.g., V. S. Berezinsky, S. V. Bulanov, V. A. Dogiel, V. L. Ginzburg, V. S. Ptuskin, 
Astrophysics of Cosmic Rays (North-Holland, Amsterdam, 1990); T. K. Gaisser, Cosmic Rays and Particle Physics, Cambridge 
University Press (Cambridge, 1998). 

59. J. R. Hoerandel, Astropart. Phys. 21, 241 (2004) [arXiv:astro-ph/0402356]. 

60. V. S. Berezinsky, S. I. Grigorieva and B. L Hnatyk, arXiv:astro-ph/0403477. 

61. A. A. Watson, arXiv:astro-ph/0408110, and references therein. 

62. from the review of particle physics, http://pdg.lbl.gov/2004/reviews/cosmicrayrpp.ps. 

63. for a recent overview see, e.g., R. Battiston, arXiv:astro-ph/0208108. 

64. W. S. Burgett and M. R. O'Malley Phys. Rev. D 67, 092002 (2003) [arXiv:hep-ph/0301001]. 

65. M. Teshima et al., Proc. 28th International Cosmic Ray Conference, Tsukuba, Japan, 1 (2003) 437, see http://www-rccn.icrr.u- 
tokyo. ac.jp/icrc2003/PROCEEDINGS/PDF/ 1 1 0.pdf. 

66. C. B. Finley and S. Westerhoff, Astropart. Phys. 21, 359 (2004) [arXiv:astro-ph/0309159]. 

67. C. B. Finley et al., Proc. 28th International Cosmic Ray Conference, Tsukuba, Japan, 1 (2003) 433, see http;//www-rccii.icrr.u- 
tokyo.ac.jp/icrc2003/PROCEEDINGS/PDF/109.pdf. 

68. K. Greisen, Phys. Rev. Lett. 16, 748 (1966); G. T. Zatsepin and V. A. Kuzmin, JETP Lett. 4, 78 (1966) [Pisma Zh. Eksp. Teor. 
Fiz. 4, 114 (1966)]. 

69. F W. Stecker, Phys. Rev Lett. 21, 1016 (1968). 

70. D. R. Bergman, Proc. 28th International Cosmic Ray Conference, Tsukuba, Japan, 1 (2003) 397, see http://www-rccn.icrr.u- 

tokyo. ac.jp/icrc2003/PROCEEDINGS/PDF/ 1 OO.pdf. 

71. J. W. Cronin, arXiv:astro-ph/0402487. 

72. J. W. Cronin, Nucl. Phys. B (Proc. Suppl.) 28B (1992) 213; The Pierre Auger Observatory Design Report (ed. 2), March 1997; 
see also http://www.auger.org. 

73. for recent reviews see, e.g., F. Halzen and D. Hooper, Rept. Prog. Phys. 65, 1025 (2002) [arXiv:astro-ph/0204527]; 
A. B. McDonald, C. Spiering, S. Schonert, E. T. Keams and T. Kajita, Rev. Sci. Instrum. 75, 293 (2004) [arXiv:astro- 
ph/03 11343]. 

74. for recent short reviews see, e.g., H. J. Volk, arXiv:astro-ph/0401122; H. J. Volk, arXiv:astro-ph/0312585. 

75. K. Mannheim, R. J. Protheroe and J. P Rachen, Phys. Rev. D 63, 023003 (2001) [arXiv:astro-ph/9812398]; J. P Rachen, 
R. J. Protheroe and K. Mannheim, arXiv:astro-ph/9908031. 

76. V. S. Berezinsky and A. Yu. Smirnov, Ap. Sp. Sci. 32, 461 (1975). 

77. P S. Coppi and F. A. Aharonian, Astrophys. J. 487, L9 (1997) [arXiv:astro-ph/9610176]. 

78. R Sreekumar et al, Astrophys. J. 494, 523 (1998) [astro-ph/9709257]. 

79. U. Keshet, E. Waxman and A. Loeb, JCAP 0404, 006 (2004) [arXiv:astro-ph/0306442]. 

80. D. F. Torres and L. A. Anchordoqui, Rept. Prog. Phys. 67, 1663 (2004) [arXiv:astro-ph/0402371]. 

81. for a review see, e.g., P Bhattacharjee and G. Sigl, Phys. Rept. 327, 109 (2000) [arXiv:astro-ph/981 101 1]. 

82. V. Berezinsky M. Kachelriess and A. Vilenkin, Phys. Rev. Lett. 79, 4302 (1997) [arXiv:astro-ph/9708217]. 

83. M. Kachelriess and D. V. Semikoz, Phys. Lett. B 577, 1 (2003) [arXiv:astro-ph/0306282]; H. B. Kim and P Tinyakov, 
Astropart. Phys. 21, 535 (2004) [arXiv:astro-ph/0306413]. 

84. T. J. Weiler, Phys. Rev Lett. 49, 234 (1982). Astrophys. J. 285, 495 (1984). Astropart. Phys. 11, 303 (1999) [hep-ph/97 10431]. 
D. Fargion, B. Mele and A. Sails, rays," Astrophys. J. 517, 725 (1999) [astro-ph/97 10029]. S. Yoshida, G. Sigl and S. j. Lee, 



Lectures on Astroparticle Physics 



February 2, 2008 



40 



Phys. Rev. Lett. 81, 5505 (1998) [hep-ph/9808324]. 

85. Z. Fodor, S. D. Katz and A. Ringwald, Phys. Rev. Lett. 88, 171101 (2002) [hep-ph/0105064]; hep-ph/0105336; JHEP 0206, 
046 (2002) [arXiv:hep-ph/0203198]; A. Ringwald, hep-ph/01 1 1 1 12. 

86. O. E. Kalashev, V. A. Kuzmin, D. V. Semikoz and G. Sigl, Phys. Rev D 65, 103003 (2002) [arXiv:hep-ph/01 12351]. 

87. D. S. Gorbunov, P G. Tinyakov and S. V. Troitsky, Astropart. Phys. 18, 463 (2003) [arXiv:astro-ph/0206385]. 

88. S. Singh and C. P Ma, Phys. Rev D 67, 023506 (2003) [arXiv:astro-ph/0208419]. 

89. G. Gelmini and A. Kusenko, Phys. Rev. Lett. 84, 1378 (2000) [arXiv:hep-ph/9908276]. 

90. V. Berezinsky, M. Kachelriess and S. Ostapchenko, Phys. Rev. Lett. 89, 171802 (2002) [arXiv:hep-ph/0205218]. 

91. For general information see http://wvvw-glast.stanford.edu 

92. A. Loeb and E. Waxman, Nature 405, 156 (2000) [arXiv:astro-ph/0003447]; F Miniati, Mon. Not. Roy. Astron. Soc. 337, 199 
(2002) [arXiv:astro-ph/0203014]. 

93. D. S. Gorbunov, G. G. Raffelt and D. V. Semikoz, Phys. Rev D 64, 096005 (2001) [arXiv:hep-ph/0103175]. 

94. G. R. Farrar, Phys. Rev Lett. 76, 41 1 1 (1996) [arXiv:hep-ph/9603271]; D. J. H. Chung, G. R. Farrar and E. W. Kolb, Phys. 
Rev D 57, 4606 (1998) [arXiv:astro-ph/9707036]. 

95. I. F. Albuquerque et al. [E761 Collaboration], Phys. Rev Lett. 78, 3252 (1997) [arXiv:hep-ex/9604002]; A. Alavi-Harati etal. 
[KTeV Collaboration], Phys. Rev. Lett. 83, 2128 (1999) [arXiv:hep-ex/9903048]. 

96. R G. Tinyakov and L L Tkachev, Phys. Rev. D 69, 128301 (2004); D. S. Gorbunov, R G. Tinyakov, L L Tkachev and 

S. V. Troitsky, Astrophys. J. 577, L93 (2002) [arXiv:astro-ph/0204360]; R G. Tinyakov and L L Tkachev, Astropart. Phys. 18, 
165 (2002) [arXiv:astro-ph/01 11305]; P G. Tinyakov and L L Tkachev, JETPLett. 74, 445 (2001) [PismaZh. Eksp. Teor. Fiz. 
74, 499 (2001)] [arXiv:astro-ph/0102476]. 

97. Z. Fodor, S. D. Katz, A. Ringwald and H. Tu, Phys. Lett. B 561, 191 (2003) [arXiv:hep-ph/0303080]. 

98. G. Domokos and S. Kovesi-Domokos, Phys. Rev Lett. 82, 1366 (1999) [arXiv:hep-ph/98 12260]. 

99. N. Arkani-Hamed, S. Dimopoulos and G. R. Dvali, Phys. Lett. B 429, 263 (1998) [arXiv:hep-ph/9803315]; I. Antoniadis, 
N. Arkani-Hamed, S. Dimopoulos and G. R. Dvali, Phys. Lett. B 436, 257 (1998) [arXiv:hep-ph/9804398]; N. Arkani-Hamed, 
S. Dimopoulos and G. R. Dvali, Phys. Rev D 59, 086004 (1999) [arXiv:hep-ph/9807344]. 

100. E. J. Ahn, M. Cavagha and A. V. Olinto, Phys. Lett. B 551, 1 (2003) [arXiv:hep-th/0201042]. 

101. J. L. Feng and A. D. Shapere, Phys. Rev. Lett. 88, 021303 (2002) [arXiv:hep-ph/0109106]. 

102. S. Hannestad and G. G. Raffelt, Phys. Rev. D 67, 125008 (2003) [Erratum-ibid. D 69, 029901 (2004)] [arXiv:hep- 
ph/0304029]. 

103. M. Casse, J. Paul, G. Bertone and G. Sigl, Phys. Rev Lett. 92, 111102 (2004) [arXiv:hep-ph/0309173]. 

104. M. Kachelriess and M. Plumacher, Phys. Rev D 62, 103006 (2000) [arXiv:astro-ph/0005309]. 

105. D. A. Morris and A. Ringwald, Astropart. Phys. 2, 43 (1994) [arXiv:hep-ph/9308269]. 

106. C. Tyler, A. V. Olinto and G. Sigl, Phys. Rev D 63, 055001 (2001) [arXiv:hep-ph/0002257]. 

107. see, e.g., L. A. Anchordoqui, J. L. Feng, H. Goldberg and A. D. Shapere, Phys. Rev. D 65, 124027 (2002) [arXiv:hep- 
ph/0 112247]; L. A. Anchordoqui, J. L. Feng, H. Goldberg and A. D. Shapere, Phys. Rev. D 66, 103002 (2002) [arXiv:hep- 
ph/0207139]. 

108. A. Kusenko and T. J. Weiler, Phys. Rev. Lett. 88, 161101 (2002) [arXiv:hep-ph/0106071]. 

109. H. Sato and T. Tati, Prog. Theor. Phys. 47, 1788 (1972); D. A. Kirzhnits and V. A. Chechin, Sov. J. Nucl. Phys. 15, 585 
(1972); L. Gonzalez-Mestres, arXiv:hep-th/0208064. 

110. S. R. Coleman and S. L. Glashow, Phys. Lett. B 405, 249 (1997) [arXiv:hep-ph/9703240]; S. R. Coleman and S. L. Glashow, 

Phys. Rev D 59, 116008 (1999) [arXiv:hep-ph/9812418]. 

111. Phys. Lett. B 497, 265 (2001) [arXiv:hep-ph/0006210]. 

112. D. Colladay and V. A. Kostelecky, Phys. Rev D 58, 116002 (1998) [arXiv:hep-ph/9809521]. 

1 13. R. C. Myers and M. Pospelov, Phys. Rev Lett. 90, 21 1601 (2003) [arXiv:hep-ph/0301 124]. 

114. see, e.g., J. R. EUis, N. E. Mavromatos and D. V. Nanopoulos, Phys. Rev D 62, 084019 (2000) [arXiv:gr-qc/0006004]. 

115. S. G. NibbeUnk and M. Pospelov, arXiv:hep-ph/0404271. 

116. R. Aloisio, P Blasi, R L. Ghia and A. F Grillo, Phys. Rev. D 62, 053010 (2000) [arXiv:astro-ph/0001258]. 

117. see, e.g., R. J. Protheroe and H. Meyer, Phys. Lett. B 493, 1 (2000) [arXiv:astro-ph/0005349]; T. Jacobson, S. Liberati and 

D. Mattingly, Phys. Rev D 66, 081302 (2002) [arXiv:hep-ph/01 12207]. 

118. J. Magueijo and L. Smolin, Phys. Rev D 67, 044017 (2003) [arXiv:gr-qc/0207085]. 

1 19. T. A. Jacobson, S. Liberati, D. Mattingly and F W. Stecker, Phys. Rev Lett. 93, 021 101 (2004) [arXiv:astro-ph/0309681]. 

120. T. Jacobson, S. Liberati and D. Mattingly, arXiv:hep-ph/0407370. 

121. for a non-techrrical overview see M. Pospelov and M. Romalis, Physics Today July 2004, p. 40. 

122. for reviews see, e.g., P. P Kronberg, Rept. Prog. Phys. 57, 325 (1994); D. Grasso and H. R. Rubinstein, Phys. Rept. 348, 163 
(2001) [arXiv:astro-ph/0009061]. 

123. J. P Vallee, Fundamentals of Cosmic Physics 19 (1997) 1 ; J. L. Han and R. Wielebinski, arXiv:astro-ph/0209090. 

124. D. Ryu, H. Kang, and P L. Biermann, Astron. Astrophys. 335 (1998) 19. 

125. P Blasi, S. Buries and A. V. Olinto, Astrophys. J. 514, L79 (1999) [arXiv:astro-ph/9812487]. 

126. G. A. Medina-Tanco and T. A. Ensshn, Astropart. Phys. 16, 47 (2001) [arXiv:astro-ph/001 1454]. 

127. E. Waxman and J. Miralda-Escude, Astrophys. J. 472, L89 (1996) [arXiv:astro-ph/9607059]. 

128. G. Sigl, M. Lemoine and R Biermann, Astropart. Phys. 10, 141 (1999) [arXiv:astro-ph/9806283]. 

129. C. Isola, M. Lemoine and G. Sigl, Phys. Rev D 65, 023004 (2002) [arXiv:astro-ph/0104289]. 

130. M. Lemoine, G. Sigl and P. Biermann, arXiv:astro-ph/9903124. 



Lectures on Astroparticle Physics 



February 2, 2008 



41 



131. T. Stanev, arXiv:astro-ph/0303123; T. Stanev, D. Seckel and R. Engel, Phys. Rev. D 68, 103004 (2003) [arXiv:astro- 
ph/0108338]; T. Stanev, R. Engel, A. Mucke, R. J. Protheroe and J. P Rachen, Phys. Rev. D 62, 093005 (2000) 

[arXiv:astro-ph/0003484]. 

132. C. Isola and G. Sigl, Phys. Rev. D 66, 083002 (2002) [arXiv:astro-ph/0203273]. 

133. R Sommers, Astropart. Phys. 14, 271 (2001) [arXiv:astro-ph/0004016]. 

134. P Blasi and D. De Marco, Astropart. Phys. 20, 559 (2004) [arXiv:astro-ph/0307067]. 

135. R Blasi and A. V. Olinto, Phys. Rev. D 59, 023001 (1999) [arXiv:astro-ph/9806264]. 

136. H. Yoshiguchi, S. Nagataki, S. Tsubaki and K. Sato, Astrophys. J. 586, 1211 (2003) [Erratum-ibid. 601, 592 (2004)] 
[arXiv:astro-ph/0210132]; H. Yoshiguchi, S. Nagataki and K. Sato, Astrophys. J. 592, 311 (2003) [arXiv:astro-ph/0302508]; 
H. Yoshiguchi, S. Nagataki and K. Sato, cosmic rays above 10**19-eV with modifications by the galactic magnetic Astrophys. 
J. 596, 1044 (2003) [arXiv:astro-ph/0307038]. 

137. R. Aloisio and V. Berezinsky, arXiv:astro-ph/0403095. 

138. G. Medina Tanco, "Cosmic magnetic fields from the perspective of ultra-high-energy cosmic rays propagation", 
Lect. Notes Phys. 576 (2001) 155. 

139. G. Sigl, F. Miniati and T. A. Ensslin, Phys. Rev. D 68, 043002 (2003) [arXiv:astro-ph/0302388]; G. Sigl, R Miniati and 
T. A. Ensslin, arXiv:astro-ph/0309695; G. Sigl, P. Miniati and T. A. Ensslin, arXiv:astro-ph/0401084; to appear in Phys. Rev. 
D. 

140. K. Dolag, D. Grasso, V. Springel and I. Tkachev, JETP. Lett. 79, 583 (2004) [aTXiv:astro-ph/0310902]. 

141. R. M. Kulsrud, R. Cen, J. R Ostriker and D. Ryu, Astrophys. J. 480 (1997) 481. 

142. For general information see http://www.euso-mission.org. 

143. E. Armengaud, F. Miniati, G. Sigl et al., in preparation. 

144. G. Bertone, C. Isola, M. Lemoine and G. Sigl, Phys. Rev. D 66, 103003 (2002) [arXiv:astro-ph/0209192]. 

145. see, e.g., L. N. Epele and E. Roulet, JHEP 9810, 009 (1998) [arXiv:astro-ph/9808104]. 

146. G. Sigl, arXiv:astro-ph/0405549; to appear in JCAP 

147. T. Yamamoto, K. Mase, M. Takeda, N. Sakaki and M. Teshima, Astropart. Phys. 20, 405 (2004) [arXiv:astro-ph/03 12275]. 



Lectures on Astroparticle Physics 



February 2, 2008 



42 



